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Double bubbles minimize 

By Joel Hass and Roger Schlafly* 
Abstract 

The classical isoperimetric inequality in M 3 states that the surface of small- 
est area enclosing a given volume is a sphere. We show that the least area 
surface enclosing two equal volumes is a double bubble, a surface made of two 
pieces of round spheres separated by a flat disk, meeting along a single circle 
at an angle of 120°. 

1. Introduction 

In this paper we find the unique surface of smallest area enclosing two 
equal volumes. The surface is called a double bubble, and is made of two pieces 
of round spheres separated by a disk, meeting along a single circle at an angle 
of 120°. This is the form assumed by two equally sized soap bubbles which are 
brought together until their boundaries conglomerate to form a common wall. 
See Figure 1, due to J. Sullivan. 

Isoperimetric problems, which study maximizing the size of an enclosed 
bounded region whose boundary is of fixed size, are among the oldest problems 
in mathematics. For a broad discussion of isoperimetric problems see Osserman 
[33]. 

The two volume isoperimetric problem in M 3 was considered by the Belgian 
physicist J. Plateau [35], and appears in C.V. Boys' famous book on soap 
bubbles. As Boys wrote, 

"When however the bubble is not single, say two have been blown 
in real contact with one another, again the bubbles must together 
take such a form that the total surface of the two spherical segments 
and of the part common to both, which I shall call the interface, is 
the smallest possible surface which will contain the two volumes of 
air and keep them separate." 



"The first author was partially supported by the National Science Foundation. 
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Figure 1. A double bubble enclosing two equal volumes. 

We now outline briefly the recent history of the problem, with more details 
given later. F. Almgren proved the existence of bubble clusters enclosing a 
given collection of volumes in three-space and minimizing area among a general 
class of surfaces. J. Taylor established a regularity result for these bubble 
clusters, showing that they satisfy regularity properties physically observed 
by Plateau. B. White suggested an argument for showing that a two-region 
minimizer must be a surface of revolution, and J. Foisy and M. Hutchings 
filled in the details of this argument. At this point it was known that the 
minimizer consisted of pieces of Delaunay surfaces meeting at 120°. Hutchings 
went on to eliminate many remaining configurations, leaving just the standard 
double bubble and an additional family known as torus bubbles as the potential 
minimizers. In this paper we closely examine the family of torus bubbles and 
show they cannot be minimizers for the equal volume double bubble problem. 
Each torus bubble is determined by an ODE, and we can rigorously eliminate 
the entire collection by an examination of the solutions of these ODEs using a 
computational technique called interval arithmetic, described in Section 5. At 
the end, only the standard double bubble remains as a possible minimizer. 

Extensive interest in the problem was generated in recent years by work 
of F. Morgan. The planar case has been solved in [1] by methods special to 
two dimensions. A fundamental paper of Almgren [3] established the existence 
of solutions to a great variety of geometric minimization problems, including 
multiple component isoperimetric problems in IR 3 . Despite the fact that exis- 
tence was established some time ago, our result is the first explicit example of 
a surface in M 3 solving a multiple region finite volume isoperimetric problem. 
Lawlor and Morgan have shown that the cone over the regular tetrahedron 
is a minimizing surface among all those with its boundary that separate the 
four faces of the tetrahedron, but it is still not known whether this surface is 
minimizing without the separation condition. 
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Such multiple region isoperimetric problems arise in many fields, including 
the growth and shape of biological cells [28], [40]. There were extensive stud- 
ies of such problems by physical and biological scientists in the 19 th century 
[24], [35]. Plateau established experimentally that a soap bubble cluster is a 
piecewise-smooth surface having only two types of singularities. The first type 
of singularity occurs when three smooth surfaces come together along a smooth 
triple curve at an angle of 120°. The second type of singularity occurs when 
six smooth surfaces and four triple curves converge at a point. The angles are 
then equal to those of the cone over the 1-skeleton of a regular tetrahedron. A 
mathematical proof that these types of singularities are the only ones possible 
in a minimizing bubble in R 3 was given by Taylor [39]. However there were 
no explicit minimizing bubbles known for any collection of volumes exhibiting 
either of these singularities. Thus we have found the first explicit example of 
a closed minimizing surface in M 3 known to exhibit some of the singularities 
predicted by Plateau. 

The arguments presented in this paper are a mixture of geometrical anal- 
ysis and of estimates of geometric quantities obtained by the use of numerical 
computation. We perform these calculations with strict estimates on the accu- 
racy of the computations. Since it is still somewhat unusual in a mathematical 
proof to use digital computers to do calculations involving real numbers, we 
will say a few words about the nature of this part of the argument. 

Computers are widely known in the mathematical community to have 
been used successfully for analysis of discrete and combinatorial problems, 
such as the 4-color theorem [6], [37]. Problems such as isoperimetric inequali- 
ties are qualitatively different, since they inherently involve real numbers. Real 
numbers are represented in digital computers in a floating point format which 
allows exact description of only a finite number of rationals. Simple calcula- 
tions, like division, or even addition, lead to unrepresentable numbers, and so 
the computer must round off. For purposes of mathematical proof, the size of 
the round-off must be accurately tracked throughout the calculation. 

Methods for strictly estimating solutions of differential equations exist, 
but are not yet widely used in the mathematical community. The vast ma- 
jority of numerical work is for approximation and simulation that does not 
meet the standards of a mathematical proof. However there have been some 
important results achieved through rigorous use of floating point numerical 
methods to achieve traditional mathematical proofs, notably Lanford's work 
on the Feigenbaum Conjectures [27] and Fefferman and de la Llave's work 
on the stability of matter [12], [15] and the work of McKay and Percival on 
nonexistence of invariant tori [34]. See also recent work of Hales on the Ke- 
pler Conjecture [19] and Gabai, Meyerhoff and Thurston [17] on hyperbolic 
3-manifolds. Computer calculations are essential to our proof that equal vol- 
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ume double bubbles minimize area. There are too many calculations to be 
done by hand. 

The proof parameterizes the space of possible solutions by a two-dimensional 
rectangle, one dimension corresponding to an angle and the other to a mean 
curvature. This rectangle is subdivided into 15,016 smaller rectangles which 
are investigated by calculations involving a total of 51,256 numerical integrals. 
Every calculation is done with strict error bounds, and all results are precise 
mathematical statements. All operations conform to the IEEE 754 standard for 
computer arithmetic, a widely adopted standardized method for implementing 
real number (floating point) computations on computers [5]. Our methods in- 
dicate that numerical techniques are likely to play an important role in future 
geometrical arguments. 

The main result is the following, announced in [20]: 

Theorem 1. The unique surface of least area enclosing two equal vol- 
umes in R 3 is a double bubble. 

Corollary 2. For any surface in M 3 enclosing two regions, each of 
volume v, the area a satisfies 

a 3 > 243tto 2 , 

with equality if and only if the surface is a double bubble enclosing two regions 
of volume v. 

The first issue in establishing Theorem 1 is to find an appropriate category 
of surface in which to minimize area. 

It suffices in this paper to consider piecewise smooth two-dimensional 
surfaces. It is sometimes useful to consider a much more general notion of 
surface, such as the (F,e,5) sets described in [3]. Our arguments actually 
imply that the double bubble minimizes in this larger class; see Theorem 2.2. 
Consideration of such a larger class of surfaces is needed primarily in the 
establishment of the existence and regularity of a minimizer, carried out in 
[3] and [39], and we will not need to be overly concerned with it in this paper. 

Define a piecewise-smooth curve to be an embedded finite union of smooth 
curves, with any two either disjoint or having intersection contained in their 
endpoints. Define a piecewise-smooth surface to be an embedded finite union 
of smooth surfaces with piecewise-smooth boundary curves, with any two sur- 
faces either disjoint or intersecting along piecewise-smooth curves contained in 
their boundaries. The set of points which are not in the interior of a smooth 
subsurface of a piecewise-smooth surface is called the singular set. 

Define a bubble to be a piecewise-smooth surface satisfying: 

1. Each two dimensional surface has constant mean curvature. 
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2. The singular set is of the type described by Plateau. It consists of smooth 
triple curves along which three smooth surfaces come together at an angle 
of 120° and isolated vertices where six smooth surfaces and four triple 
curves converge at a point. The angles at the point are equal to those of 
the cone over the 1-skeleton of a regular tetrahedron. 

3. The mean curvatures around an edge in the singular set where three 
surfaces have common boundary sum up to zero. 

The above conditions are necessary for no local perturbation of a piecewise- 
smooth surface to decrease the area while preserving the volume in each of its 
complementary regions. 

A bubble enclosing regions of prescribed volumes is called a minimizing 
bubble if it minimizes area among all bubbles enclosing the same volumes. The 
regions are not necessarily connected. 

Given positive constants v\ and V2, let a(v±, V2) denote the infimum of the 
area of piecewise-smooth surfaces enclosing two regions i?i and R 2 in M 3 which 
are closed bounded sets with disjoint interiors such that volume(i?i) = v\ and 
volume(i?2) = v 2- We will refer to the interior of the complement of R\ U R 2 
as the exterior region Rq. Note that i?o may not be connected, in which case 
the bubble encloses some compact "empty regions" . 

A double bubble enclosing volumes v\,V2 is a surface made of three pieces of 
round spheres, meeting along a single circle at an angle of 120°, and enclosing 
two connected regions having volumes v\ and v 2 . We consider the plane to 
be a sphere of infinite radius in this setting, allowing the interface of a double 
bubble to be a fiat disk, as occurs when v\ and v 2 are equal. See Figure 2. 
In this and all other figures, we assume that the axis of rotational symmetry 
is the x-axis. Furthermore, we will always take the generating curves for a 
surface of revolution to be in the upper half of the xy-plane. 

It has been conjectured since the work of Plateau that double bubbles give 
the most efficient shape for enclosing two given volumes. 

Conjecture 3. The double bubble enclosing volumes v\,v 2 is the unique 
surface having area equal to a(vi,v 2 ). 

Theorem 1 solves this conjecture in the case that v\ = v 2 - 
Deep results of Almgren and Taylor, summarized in Section 2, imply that 
for any two positive numbers v±, v 2 there exists a minimizing bubble S(vi,v 2 ) 
in M 3 which encloses volumes v\ and v 2 . Arguments based on ideas of B. White 
and F. Morgan, and developed in [16], [22], [32], show that any solution must 
be a surface of revolution. We will need to refer to this proof, so we present a 
simple version for the case of two regions in M 3 in Theorem 2.6. The lack of 
such an argument for isoperimetric problems involving three or more regions 
in M 3 makes those problems more formidable. 
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While the reduction to a surface of revolution gives an enormous simpli- 
fication in the scale of the problem, F. Morgan has pointed out several major 
topological obstacles to solving the double bubble conjecture. The first is that 
the regions Ri and R2 bounded by the minimizing bubble S(v±, 1*2) may not 
be connected, as illustrated in Figure 3. 

A second problem is that S(vi,V2) may enclose bounded regions in R 3 
which do not form part of either R\ or R2, as in Figure 4. These regions are 
called empty regions by Morgan. 




Figure 2. Equal and unequal volume double bubbles in cross-section. 
The axis of revolution is the x-axis. 




Figure 3. Cross-section of a bubble with a nonconnected region. 
The nonconnected region is shaded. 




Figure 4. Cross-section of a bubble with an empty torus region. The 
shaded region generates a solid torus when revolved around the x-axis. 
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A third problem is that S(vi,V2) may enclose nonsimply connected re- 
gions. Since the surface is a surface of revolution, these regions are homeomor- 
phic to solid tori, and we call them torus components. 

There are therefore numerous possible configurations which a minimizing 
bubble might take. A recent breakthrough due to Hutchings [22] has given 
restrictions on the type of surfaces that can arise in the double bubble problem, 
eliminating many possibilities. Some remaining possibilities are depicted in 
Figure 5. For the case of equal volumes, Hutchings showed that there were 
further constraints. Each region must be connected, leaving only two possible 
configurations, the double bubble and an additional class of possibilities, called 
torus bubbles, whose properties will be discussed in Section 4. 




Figure 5. Cross-sections of some possible bubble configurations. 

In Section 5 we describe the algorithm used in a series of computations 
which show that torus bubbles are not minimizers for the two equal volume 
isoperimetric problem in R 3 . 

Techniques extending those developed in this paper also prove Conjec- 
ture 3 for other volume ratios, but at this time do not seem to suffice to cover 
all cases. 



2. Existence and regularity 

Almgren showed in [3] that there exists an area minimizing surface S(vi,V2) 
in M 3 among the set of surfaces enclosing volumes v\ , V2 ■ Here surface refers to 
a generalized notion defined using the methods of geometric measure theory, 
specifically what Almgren calls (F, e, 5) sets. For our purposes it suffices that 
this class includes the piecewise-smooth surfaces. 
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Theorem 2.1 ([3]). An area minimizing surface S(v±,V2) exists and is 
a smooth surface almost everywhere. 

The nature of the singularities of S(vi, V2) was established by Taylor. 

Theorem 2.2 ([39]). S (1)1,1)2) is a piecewise- smooth surface. Its singu- 
larities consist of smooth triple curves along which three smooth surfaces come 
together at an angle of 120° and isolated points where four triple curves and 
six pieces of surface converge. At these isolated points the asymptotic cone is 
the cone over the 1- skeleton of a regular tetrahedron. 

Thus Taylor's work established that S(v±, V2), which Almgren's theorem 
established as a minimizer among a generalized class of surfaces called (F, e, 5) 
surfaces, is a bubble, in our terminology. It follows that if we can establish the 
minimality of the double bubble among our class of bubbles then we will also 
show it minimizes among the more general class of surfaces considered by Alm- 
gren. We now establish some properties of minimizing bubbles. Lemma 2.3, 
Corollary 2.4 and Lemma 2.5, are standard results in variational geometry. 
We include brief proofs for completeness. Lemma 2.3 states that the sum of 
the (oriented) mean curvatures of all points crossed by a path starting and 
ending in the same region is zero. It implies that the sum of the (oriented) 
mean curvatures of all points crossed by any path depend only on its starting 
and ending points. These results hold even if the path crosses through many 
surfaces, and if it starts and ends in different components of the same region. 

Lemma 2.3. Let 7 be an oriented curve in M 3 intersecting a minimizing 
bubble B transversely at regular points such that the initial and final points of 
7 lie in the interior of the same region. Then the sum of the mean curvatures 
of all the points of 7 fl B, oriented by 7, is zero. 

Proof. Perturb the curve slightly so that each of its intersections with B 
becomes perpendicular. Consider a deformation of B which pushes points in 
R 3 near 7 a uniform distance along the curve. To first order this preserves the 
volume of each region. The derivative of the area, to first order, is given by the 
sum of the mean curvatures over the points of 7 n B. If this sum is nonzero, a 
deformation can be defined which decreases area while preserving the volume 
of each region. □ 

A special case of the above lemma occurs when 7 is a simple closed curve 
encircling a triple curve of the bubble. The lemma then implies that the sum 
of the mean curvatures around the triple curve adds up to zero. This local 
minimization condition is built into our definition of a bubble. More generally, 
the lemma implies that any two surfaces separating the same pair of regions 
have the same mean curvature. 
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Corollary 2.4. // S± and S2 are surfaces in a bubble, each of which 
separate regions R\ and R2, then their mean curvatures are equal. 

Proof. Consider an arc a which starts in R\ passes through S\ into R2, 
then transversely through the bubble and on through S2 back into R\. A 
subarc (3 of a starts and ends in R 2 . Applying Lemma 2.3 to both a and (3 
implies the corollary. □ 

Lemma 2.5. A minimizing bubble B is connected. 

Proof. If not, we can move a component by an isometry of M 3 until it 
touches a distinct component. The resulting singularity violates those allowed 
in Theorem 2.2. □ 

The following result was first observed by B. White and F. Morgan, and 
first written down by Foisy in [16]. Since we will need to refer to the proof, 
we present an argument for the case of two regions in R 3 . More detailed and 
general arguments can be found in the work of Hutchings [22] and Morgan [32]. 

Theorem 2.6. S{v\,V2) is a surface of revolution. 

Proof. Since S(vi,v 2 ) is a minimizing bubble, each of its faces separates 
a pair of distinct regions - either separating R\ from R2 or one of them from 
Rq (which is not connected if there are empty regions.) 

Given any unit vector Z in R 3 , there is a plane Pz in R 3 perpendicular to Z 
which bisects the total volume v\ + V2 enclosed by S(v\,v 2 ). If Qz is a parallel 
plane which also bisects the total volume, then Pz and Qz split the regions R\ 
and R 2 in the same proportions. (In fact Pz must coincide with Qz, but we do 
not need to show this. It is also irrelevant to this argument whether S(vi,v 2 ) 
contains empty regions.) Let f(Z) denote the proportion of the volume of 
R\ on the side of Pz to which Z points. The function f(Z) is well defined, 
continuous on the unit 2-sphere, < f(Z) < 1 and f{—Z) = 1 — f(Z). Along 
any great circle of directions on the unit 2-sphere, the Intermediate Value 
Theorem implies that there are at least two points where f(Z) = 1/2 and the 
plane Pz bisects the volume of both regions. We fix Z to be such a vector. 

Consider the intersection of S{v\,v 2 ) with each of the two half-spaces 
determined by Pz- The intersection of a face of S(vi,v 2 ) with the plane Pz 
is given by a graph in Pz with smooth edges and isolated vertices, or else the 
face is completely contained in Pz, since the face is a smooth constant mean 
curvature surface [18]. Reflection of the smaller area piece of S(vi,v 2 ) lying in 
one of these half-spaces gives a new surface Si(vi,v 2 ) which encloses regions 
of the same volumes. Since <Si(i>i,i>2) cannot have less area, it must have the 
same area as S(v\,v 2 ). S\(vi,V2) has the property that reflection through Pz 
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preserves the new regions, which we continue to call Ro, R\ and i?2- If the 
plane Pz contained a face of S(v\,V2) then this face would now separate two 
components of the same region in Si(v±,V2), and area could be reduced by 
removing that face from S±(vi, V2), preserving the volume of each region. So 
Pz cannot contain a face of S(y\, ^2). 

We repeat the above argument along a great circle of the unit 2-sphere 
consisting of the unit vectors which are perpendicular to Z, to find a plane 
P\Y perpendicular to Pz which bisects the volume of the regions Ri and R2 
in Si(v±, U2). This process constructs a surface S^C^i, ^2), still having the 
same area as S(v±,V2), enclosing regions of the same volumes, and for which 
reflection through each of Pz and Pw preserves R\ and i?2- 

Since composing reflections through two perpendicular planes gives a ro- 
tation of angle 180°, it follows that rotation of angle 180° about L = Pz l~l Pw 
preserves the regions i?i and R2. Now consider any plane Q containing L. 
Note that half of the volume of R\ and i?2 lies on each side of Q. As before 
the plane Q cannot contain a face of S^^i, ^2)- If the intersection of Q with 
£2(^1,^2) is not perpendicular to Q, then replacing half of £'2(^1,^2) with its 
reflected image through Q gives a new minimizing surface which is not smooth, 
and has singularities not allowed by Theorem 2.2. Thus ^(fi, ^2) is perpen- 
dicular to each plane through L and is therefore a surface of revolution around 
L. <Si(t>i,t>2) and S2(vi,V2) may not a priori be identical, but they coincide 
on a half-space of M 3 . Moreover this half-space was chosen arbitrarily, so that 
<Si(^ij ^2) is cut by Pw into two pieces, each half of some surface of revolution. 
The axis of each of these surfaces of revolution is L = Pz n Pw, so they co- 
incide and Si(v±,V2) itself is a surface of revolution. Similar reasoning shows 
that S(vi,V2) is a surface of revolution, and coincides with S^^i,^) every- 
where. If we make Si by reflecting the other half of S across Pz, then the same 
argument shows that the other half of S is also a surface of revolution about 
a line L' in Pz- The Almgren and Taylor regularity results, and in particular 
the unique continuation property of constant mean curvature surfaces, imply 
that V = L, so S itself is a surface of revolution about L. □ 

We next summarize some key results obtained by Hutchings in [22]. Recall 
that 5(^1,^2) is a minimizing bubble that separates M 3 into regions R\ of 
volume v\ and R2 of volume V2 and that a(vi,V2) is the area of S{v\,vz). 

Theorem 2.7. The function a{v\,V2) is strictly concave on [0, +00) x 
[0,+oo). 

Corollary 2.8. Ifv[ > v\ then a(v[,V2) > a(vi,v 2 ). 



Proof. As v\ — ► +00, 0(^1,^2) > a(v\) — ► +00. Concavity of the function 
o implies the corollary. □ 
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Corollary 2.9. S(vi,V2) has no empty regions. 

Proof. If there is an empty region, join it it to one of the other regions by 
removing a face on its boundary and apply the previous corollary. □ 

Deeper connectedness results also follow from concavity. In particular, 
Hutchings deduced that each of Ri and R2 is connected if the volumes are 
equal. 

Theorem 2.10. S(vi,v±) encloses exactly two connected components. 
Proof. See Hutchings [22, Th. 4.2]. □ 



3. Delaunay surfaces 

In this section, we summarize the classification theory of constant mean 
curvature surfaces of revolution, the Delaunay surfaces, and present some prop- 
erties of these surfaces that will be used in our study of bubbles. 

The mean curvature of a surface in K 3 is the trace of the second fundamen- 
tal form of the surface. Computing this requires choosing a unit normal vector 
field to the surface. The mean curvature vector field, formed by scaling a unit 
normal vector field by the mean curvature, gives the direction of a variation 
which decreases the area of a surface as quickly as possible [38]. 

For a surface of revolution we have simple formulas for the mean curvature. 
Consider a surface of revolution about the x-axis in M 3 , with generating curve 
contained in the upper half-plane. If the generating curve is a graph, y = y(x), 
then the mean curvature h, equal to the sum of the principal curvatures, is 
given by the formulas: 

Km (1 + y2)3/2 > 

k - 1 

yv l + y 

h — k m -\- kp, 

where k m is the curvature at (x, y) of the generating curve in the xy-plane 
(sometimes called the meridian, or profile curve) and k p is the normal curva- 
ture of the parallel curve. The value of k p is equal to the reciprocal of the 
distance to the x-axis along the perpendicular to the generating curve [38]. 
Note that the mean curvature of a unit sphere equals two. For surfaces of 
revolution generated by curves which are not graphs, we need to specify an 
orientation on the generating curve to fix the sign of the mean curvature. We 
use the convention that the sign of the mean curvature is given by the formulas 
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above if the curve is oriented left-to-right, or equivalently if the graph's tan- 
gent vector has a positive x-component. If its tangent vector has a negative 
x-component, the signs in the formulas for k m and k p should be reversed in the 
above formulas. For curves with vertical tangent vectors the above formulas do 
not apply directly, but in our applications vertical tangents are isolated points, 
and the limits of the above formulas will give k m and k p at vertical tangencies. 

Surfaces of revolution having constant mean curvature were first studied 
by Euler, and classified by Delaunay [13]. The term Delaunay surface is used 
to refer to these surfaces. 

Theorem 3.1. The surfaces of revolution of constant mean curvature 
are the plane, sphere, catenoid, cylinder, unduloid, and nodoid. The generating 
curves of a nodoid and an unduloid, called the nodary and the undulary, are 
periodic along the x-axis, and have exactly one local minimum and one local 
maximum in each period. The undulary is a graph over the x-axis. The nodary 
has one local maximum, one local minimum and two vertical tangencies in each 
period. 

Proof. Expositions of the classification of Delaunay surfaces can be found 



The nodoid and unduloid are not as well known as the other surfaces. An 
undulary can be obtained by rolling an ellipse along the x-axis and tracing the 
path taken by one of the foci. A nodary can be obtained by rolling a hyperbola 
[14]. See Figure 6. 



Figure 6. Portions of a nodary and an undulary, generating curves of a 
nodoid and an unduloid. 

With their generating curves oriented left to right, the mean curvature 
of each of the Delaunay surfaces is nonnegative, with the exception of the 
nodoid. The sign of the mean curvature of a nodoid depends upon the choice 
of orientation for its generating nodary, and can be positive or negative. 

If a surface of revolution minimizes area among all surfaces of revolution 
surrounding a given volume, then its generating curve satisfies an associated 
Euler-Lagrange equation. This second order ODE implies that the mean cur- 
vature is constant [14]. 



in [14] and [23]. 



□ 
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The mean curvature of the surface of revolution generated by a C 2 graph 
y(x) is given by: 

W h = k m + k p = JL3/0 + 



(1 + 2/ 2 ) 3 / 2 Wi+F ' 

We can integrate this ODE once, after multiplying through by 2yy. 



h • 1 



2hyy + 

df(x) 
dx 



(l+y2)3/2 y ^l^2 

Zyyy %y 



(i + y2)3/2 yrq:^ 



where /(x) is the function 
(2) f(x) = hy 



2 % 



Therefore / is constant along the graph. The fact that the equation for con- 
stant mean curvature surfaces of revolution has a first integral was known to 
Plateau [35, pp. 138-139], who in turn references work of Beer. Korevaar, 
Kusner and Solomon called / the force of a constant mean curvature surface, 
and pointed out that it has a physical interpretation as the net force exerted 
by a soap film on a plane cutting off an end of the surface. They also showed 
that it could be defined for more general constant mean curvature surfaces 
[25], where it serves as a useful analytical tool. 

A convenient way of expressing the force of a constant mean curvature 
surface of revolution which holds even when the generating curve is not a 
graph is given by the formula 

(3) f = hy 2 -2ycosa . 

Here a denotes the angle between the positive x-axis and the generating curve 
of the surface of revolution, which is an oriented curve. This makes sense for 
all values of a. In our applications, we alway orient the generating curve of a 
Delaunay surface, and thus fix an unambiguous sign for the mean curvature. 

Lemma 3.2. Given h, xq, yo > and yo, there exists a unique Delaunay 
curve 5 through the point (xo,yo) having slope yo and generating a Delaunay 
surface D having mean curvature h. The slope yo can be infinite. The Delaunay 
surface is a sphere S if and only if 



(4) 
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If h > hs, then D is a nodoid and 5 lies strictly below the circle a gen- 
erating S in a neighborhood of (xo>2/o)- An arc of 5 which is decreasing as it 
leaves (xo,yo) remains beneath a as long as it remains a graph over the x-axis. 

If h < hs, 5 lies strictly above a in a neighborhood of (xo,yo)- If < h < 
hs then D is an unduloid or cylinder, ifh<0 then D is a nodoid, and if h = 
then D is a catenoid. 

Proof. Solving equation 1 for y gives a second order ODE satisfied by the 
generating curves of Delaunay surfaces at points where y is finite: 



Given a constant h and a choice of initial conditions y and y this equation 
has a unique solution from standard existence and uniqueness of solutions to 
ODEs. These solutions can also be obtained by rolling an appropriate ellipse, 
hyperbola or parabola [14]. The latter approach also shows that we can find a 
nodoid passing through y with any mean curvature h ^ in the case where y 
is infinite. 

The mean curvature of the unique sphere centered on the x-axis and having 
slope yo at (xo,yo) is given by hs- Assume now that a is centered at (0,0) on 
the xy-plane. 

If h > hs then Equation 5 implies that 5 lies strictly below a in a neigh- 
borhood of (xo,yo). If D is an unduloid then 5 is an undulary whose graph 
crosses a again in at least two more points (x±,yi) and (£2,2/2)- We can as- 
sume, after performing a left/right reflection if necessary, that xo > and 
x\ > xo- Denote by 7 the subcurve of 5 starting at (xo,yo) and running to 
(xi, yi). Let jt = 7 + (t, 0) be a horizontal translate of 7 by a distance t and 
let T = sup{t : 74 n S* / 0}. Then jt is tangent to S at a point P and lies 
above a in a neighborhood of P. This violates Equation 5, so D cannot be an 
unduloid. D cannot be a cylinder or a catenoid, so it can only be a nodoid. 
If a subarc 7 of 5 which is decreasing as it leaves (xo,yo) intersects a while 
8 remains a graph then an identical argument shows that we can horizontally 
translate 7 until it is tangent to a while lying above it, which is impossible by 
equation 5. (In fact it will remain beneath a for somewhat longer than this, 
but this will not be relevant for us.) 

If < h < hs then equation 5 implies that S lies strictly above a in a 
neighborhood of (xo, yo)- If -D is a nodoid, consider the subarc (5 of 5 decreasing 
from (xo,2/o) until it reaches a vertical tang y v ). The arc j3 is a graph 

which may or may not cross a at an additional point. In either case we can 
translate [3 horizontally to the left until the last time at which it intersects a. 
The final intersection must be at an interior point of (5, so after translation 



>-,2 



(5) 



y = -fc(l+y2)3/2 + 



V 
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(5 becomes tangent to a at an interior point while lying underneath a, again 
impossible by equation 5. Thus the surface must be an unduloid or cylinder. 
If h < then the surface must be a nodoid. □ 

Nodoids and unduloids can be distinguished by the sign of the product 

fh. 

Corollary 3.3. For a nodoid, fh > 0. For an unduloid and a cylinder 
fh < 0. For a sphere, a catenoid and a plane, fh = 0. 

Proof. For a sphere of radius r, we can compute / at a maximum point 
where 

f = hr 2 - 2r. 

Since h = 2/r, we have / = 0. We also compute fh for an unduloid or cylinder 
at a local maximum. By Lemma 3.2 we know that < h < h$ = 2/r. Then 

/ = hr 2 - 2r < (2/r)r 2 - 2r = 0. 

So h > and / < implying that hf < 0. 

We compute fh for a nodoid at a vertical tangency, where y = +oo and 
f = hy 2 . Then fh = h 2 y 2 > 0. □ 

The next proposition deduces some useful properties of Delaunay surfaces. 

Proposition 3.4. Let D be a Delaunay surface with generating curve 5, 
mean curvature h, maximum y-value yu and minimum y-value y m . Then 

1. If 5 is a nodary or undulary, the arc-length of one period of 5 is 2?r/h. 

2. If S is an undulary, the period P of 5 satisfies 

2(yM + y m ) < P < Tr{y M + y m )- 

3. If 5 is a nodary, 

h = 2/(y M - y m )- 

If 5 is an undulary, 

h = 2/(y M +y m )- 

4. If 5 is a nodary, 5 has nonzero curvature. 

Proof. 

1. In [23] the equation for the generating curve of a Delaunay surface is 
given in terms of arclength, and it is shown that nodaries and undularies 
are periodic, with the arc length of one period given by 2ir/h. 
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2. The period of an undulary is equal to the perimeter of the ellipse that 
is rolled along the x-axis to generate it [14]. The inequalities follow 
by comparing twice the length of the major axis, 2(dm + Um), with the 
perimeter of an enclosing circle whose diameter is the length of the major 
axis, ir(y M + y m )- 

3. If 5 is a nodary, we calculate the force at yu and y m . At yu we get 

/ = hyli - 2 yM- 

At y m we get 

/ = hy 2 m + 2y m , 

where the sign of the second term changes since the orientation has re- 
versed. Setting these equal we get 

2 



h 



yu - y m 

For an unduloid, the proof is identical except for one sign change. 
At a point where S has zero curvature, k m = so 

h — kn 



y^l + y 2 

Equation 4 gives the mean curvature of a sphere with slope y at height 
y as: 

hs = A -2 > - 

yyi + y 

Where y is defined, this is exactly twice the value of h calculated above. 
For a nodary, h > hs or h < by Lemma 3.2, so we cannot have hs = 2h 
and there are no points of zero curvature on a nodary. At points where 
the nodary is vertical, h = k m since k p = 0. Thus the curvature is 
nonzero at vertical tangencies except possibly when h = 0, which would 
imply that the Delaunay surface is a plane or a catenoid, rather than a 
nodoid as assumed. 

□ 

Important stability formulas for constant mean curvature surfaces were 
developed by J. L. Barbosa and M. do Carmo [8]. Stability refers to the 
behavior of a surface when a compact subsurface is deformed by a variation, 
while holding its boundary fixed. A subsurface of a complete constant mean 
curvature surface is called stable if there is no compactly supported normal 
variation which decreases its area while preserving the volumes on each side, 
and unstable otherwise. Unstable constant mean curvature surfaces cannot 
form part of a minimizing bubble. 
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Proposition 3.5. The smooth subsurfaces of a minimizing bubble S(vi,V2) 
are stable subsurfaces of Delaunay surfaces. 

Proof. S(vi,V2) is a piecewise-smooth surface of revolution whose smooth 
subsurfaces have constant mean curvature by Theorems 2.6 and 2.2. Any 
constant mean curvature surface of revolution in M 3 is a Delaunay surface. If 
one of these subsurfaces is unstable, then there is a variation which maintains 
the volume of each region while decreasing the area. The surface then cannot 
be part of a minimizing bubble. □ 

We next deduce a formula for the horizontal distance between two points 
on a Delaunay curve whose y-coordinates are known. We require that the 
second point be above the first, but it can be either to the left or to the right. 
We also get a formula for the volume under a Delaunay surface. 

Proposition 3.6. Let (xi,y±) and (x2,yz) be two points on a Delaunay 

dy 

curve. Suppose that y\ < y 2 and that between the points, — / 0. Then 

dx 

where t = t(y) = hy 2 — f. 

dx 

If x\ < X2 and the curve also satisfies — > between the two points then 

dy 

the volume underneath the surface of revolution generated by the curve between 
x = X\ and x = x 2 is given by: 

(7) v(x 1 ,x 2 ) = r d 

^ ! ^ ! Jyi ^(2y + t)(2y-t) 

If instead we have x 2 < x\ in the previous case then the formula gives the 
negative of the volume underneath the surface of revolution generated by the 
curve between x = x\ and x = x 2 - 



Proof. We solve equation 2 for dx/dy, using that dy/dx ^ 0, and then 
apply the change of variables formula. 

/ = % 2 - 2y 



dx 



2y 
t 



dy y/(2y + t)(2y-t) 
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Note that the sign of t = t(y) = 2y cos a is positive if and only if the curve 
is oriented to the right, i.e. if cos a > in equation 3. Hence the integrand 
gives the correctly signed measure of displacement in the x direction. 

The volume of the region enclosed by a surface of revolution generated by 
a graph y{x) between x\,X2 with x\ < X2 is given by 



Since we have already deduced a formula for dx, the volume formulas follow 



Note that if y\ > y2 then the formulas in Proposition 3.6 would have the 
sign reversed. A choice in sign was made in taking the positive square root, and 
this choice is correct if the curve is oriented so that y is increasing. However 
we need to assume that X2 > x\ only when calculating the volumes, and not 
in equation 6. 

The integrals of Proposition 3.6 are singular at a local minimum or max- 
imum, where dy/dx = 0, and we need to apply a change of variables to obtain 
a formula which holds near such points. Note that the formulas below do not 
require that yi > yi, but do require that the Delaunay curve is a graph over 
the x-axis, with X2 > x±. 

Proposition 3.7. Let (x\,yi) and (x2,y2) be two points on a Delaunay 
curve which is a graph over the x-axis with X2 > x\ and suppose there is 
exactly one critical point between the two points, with y-value labeled £/min or 
Umax depending on whether it is a local minimum or maximum. Let h, f be the 
mean curvature and force associated to the Delaunay surface generated by the 
curve. Set t = t(y(z)) = h(y(z)) 2 — f. 

If Vmin is a local minimum, let y = y(z) = y min + z 2 , z 1 = -y/y[ - y m i n 
and z 3 = ^y 2 - y min . Then 




immediately. 



□ 





and 




_• 



(9) 



is a local maximum, let y(z) = y max - z 2 , z 2 = V^max - 



?/2 and 



za = -\/ymax - m- Then 



(10) 




and 



(11) 
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Finally, if the curve is an increasing graph over the x-axis between x\ and X2 
then equation 8 continues to hold if the sign of z\ is changed. 

Proof. These formulas follow from applying a change of variables to y in 
equations 6 and 7. Suppose y m \ n is a local minimum. Divide the Delaunay 
curve into the two pieces on each side of the minimum, and apply Proposi- 
tion 3.6 and the substitution y = y m i n + z 2 to each piece. 

One of the terms in the square root in the denominator in equation 6 
factors as 

2y - t = 2y - hy 2 + / = (y - y min )(2 - hy min - hy) = z 2 (2 - hy min - hy) 

since / = Ziy^m — 2y m in- A resulting factor of z in the denominator, which 
causes the singularity, cancels with a z in the numerator from the change of 
variables. The integrands are even functions of z, and y = y m - m corresponds to 
z = 0, so the curve portions on opposite sides of the minimum can be mapped 
to positive and negative z values, and combined into one integral. 

The maximum case is similar. Near a maximum the substitution y = 
2/max — z 2 is used. Several minus signs appear in the derivation, resulting in 
the given expression. 

The final statement in the proposition is a direct application of the change 
of variables formula to equation 6. □ 

4. Torus bubbles 

We established in Lemma 3.5 that a minimizing bubble must be obtained 
by revolving a union of Delaunay curves contained in the upper half plane. 
A key case of such a bubble is the torus bubble, constructed as follows. Take 
two circular arcs of the same radius, facing each other, each with one endpoint 
and center on the x-axis, and connect the other endpoints with two different 
Delaunay curves meeting at 120 degrees. Rotating around the x-axis, we get 
a piecewise-smooth surface surrounding two components, one homeomorphic 
to a torus, which we call the torus component T, and one homeomorphic to a 
ball, which we call the ball component B. It is not immediately clear whether 
it is possible to make such a construction so that the curves meet at 120° 
angles and the mean curvatures sum to zero around each triple curve. Such 
torus bubbles do indeed exist, and we will need to show that none of them are 
minimizers. Figure 7, due to J. Sullivan, shows a torus bubble. 

The reason that torus bubbles play so central a role in our argument is 
due to another key result of Hutchings: 

Theorem 4.1. A minimizing bubble enclosing two equal volumes must 
be either a double bubble or a torus bubble. 
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Figure 7. A torus bubble in R 3 . 

Proof. Theorem 2.10 states that S(v,v) contains two connected compo- 
nents. The two possible configurations are then a consequence of Hutchings' 
structure theorem for minimizing bubbles [22, Th. 5.1]. □ 

We label the surfaces of a minimizing torus bubble as indicated in Figure 8. 
The component homeomorphic to a solid torus is denoted T and the component 
homeomorphic to a ball is denoted B. The volumes of the inner ball component 
and outer torus component are denoted by V{ and v respectively. The inner 
Delaunay surface on the boundary of T is denoted Tj and the outer Delaunay 
surface is denoted T Q . The generating curve of T Q is denoted by r Q and that 
of T{ by Tf. The angles subtended by the generating curves for Si and 52 are 
denoted by 0i and 82 ■ The mean curvatures of Tj and T Q are hi and h Q , with 
signs chosen so that so that hi is positive if the mean curvature vector points 
into B and negative if it points into T, and h is positive if the mean curvature 
vector points into T and negative if it points towards outside the bubble. 

Denote by (2:1,2/1) the initial point where the two Delaunay surfaces start, 
and (2:2,2/2) t ne point where they rejoin. To specify a sign for the mean cur- 
vature, we orient the generating curves 73 and r Q of Tj and T a so that they run 
from (xi,yi) to (2:2,1/2)- 



°f T \^2-y2) 




/ \ ^ \ 
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Figure 8. Parameters of a torus bubble. 

Torus bubbles for various ranges of 61 are indicated in Figures 18, 19 
and 20. We will show that 61 and h parameterize the space of torus bubbles 
and that no choice of 61 and h Q gives a minimizing bubble. To do so, we first 
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present some properties of minimizing torus bubbles, mainly derived from the 
work of Hutchings. 

Lemma 4.2. In a minimizing torus bubble, S\ and S2 have the same 
mean curvature. 

Lemma 4.3. In a minimizing torus bubble, h Q > and hi = 2 — h Q . 

Proof. The three mean curvatures around a triple curve sum to zero by 
Lemma 2.3 implying that h Q = 2 — hi. If h Q < then there is a deformation 
which pushes T a outward, increasing the volume of the torus component and 
decreasing its area. Corollary 2.8 then gives a contradiction. □ 

Proposition 4.4. The regions bounded by a minimizing torus bubble 
are not all preserved by reflection through some plane P perpendicular to the 
x-axis. 

Proof. If there is such a symmetry, then the bubble is invariant by reflec- 
tion through the xy-plane and by reflection through a perpendicular plane P. 
The argument given in Theorem 2.6 shows that the bubble is a surface of rev- 
olution around the line of intersection of these two planes. But it is also a 
surface of revolution around the x-axis, and therefore a surface of revolution 
around two perpendicular axis. It must then be a union of spheres, and not a 
torus bubble. □ 

Corollary 4.5. Given a minimizing equal-volume torus bubble, we can 
rescale and reflect so that 9\ < 62 and each of the two circular arcs generating 
S\ and S2 has curvature one. 

Proof. If 0i = 62, Lemma 4.4 gives a contradiction. If Q\ > O2, a left /right 
reflection interchanges the two angles. Lemma 4.2 implies that the two circular 
arcs generating S\ and S2 have the same curvature, and thus the same radius. 
We can rescale and assume without loss of generality that each has radius one 
(and thus mean curvature two). □ 

From now on we will restrict attention to torus bubbles with 6\ < 62 and 
with spherical caps of radius one. 

Proposition 4.6. // a torus bubble contains a surface T\ which is a 
subsurface of a nodoid and the nodary generating T\ contains two interior 
vertical tangencies, then T\ is unstable and not a part of a minimizing bubble. 

Proof. We will construct a volume preserving Jacobi field supported on a 
proper subset of T\, establishing instability as in [8, Prop. 2.24]. 
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We start with the constant unit vector field v pointing in the direction of 
the positive y-axis. The vector field v is tangent to T\ along two circles which 
are generated by rotating the interior vertical tangencies of the nodary. Let T{ 
denote the compact portion of T\ between these two circles. Define the vector 
field W = (v,N)N. where N is the normal vector field along T[. and extend 
W to be the zero vector field elsewhere on T\. The integral of the function 
/ = (v, N) along T\ is zero, since the involution (x, y) — > (—x, —y) of T\ takes 
/ to — /. So the normal variation defined by / preserves volume in the sense 
of [8, Prop. 2.24]. Moreover [8, Prop. 2.24] shows that f'N is a Jacobi field, 
since it is constructed by taking the normal part of a constant vector field, and 
since the nodary extends beyond the two vertical tangencies, the corresponding 
Jacobi field is supported on a proper subset of T\. The Morse Index Theorem 
implies that the surface T\ is unstable and that there is a smaller area surface 
separating the same volumes. □ 

The next proposition establishes that minimizing torus bubbles do not 
contain Delaunay curves which are longer than a full period. 

Proposition 4.7. IJti orr Q contains a full period, then the torus bubble 
is not a minimizer. 

Proof. Suppose to the contrary that a minimizing torus bubble contains 
a subsurface generated by a curve which contains a full period. If a boundary 
surface of the torus component is a subsurface of a nodoid with width more 
than a period, then the associated nodary contains two vertical tangencies, 
which implies instability by Proposition 4.6. Similarly, if a boundary surface 
of the torus component is a subsurface of a cylinder of radius r and its width 
is greater than 2nr, the period of a cylinder, than the surface is shown to be 
unstable in [8]. 

If a boundary surface of the torus component is a subsurface U of an 
unduloid then the lemma is less obvious. Note that we cannot necessarily 
construct a volume preserving Jacobi field from the vector field d/dx, since the 
Delaunay curve may not contain two maxima or two minima. We establish the 
lemma with a reflection argument. Assume first that the unduloid containing 
a full period is the inner surface Tj. 

An arc u of an undulary which is longer than a period contains three 
distinct points with the same y value, (x±, y±), (x2, yi), (£3, yi), with x\ < 
X2 < x$, and with the slopes at (xi,y±), (xs,yi) each being nonzero and equal 
to the negative of the slope at (x2, yi)- Reflect the part of u between the planes 
x = x\ and x = X3 through the plane x = {x\ + x%)/2 to get a new surface of 
revolution U' generated by u', the reflection of u. Reflection preserves both the 
area and the volume under U between the planes x = x± and x = X3. Suppose 
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first that v! does not intersect t q . Then the volume of the torus bubble is 
also preserved. The reflection creates a surface which is not smooth along the 
planes x = x\ and x = X3 since reflection changes the sign of the slope of the 
reflected surface along these planes, and this slope was not zero. Thus we have 
found a solution to the minimization problem which contains a surface which 
contradicts Theorem 2.2. 

Possibly v! does intersect the outer curve r . See Figure 9. In that case the 
total volume enclosed by the union of U' and T Q between x = x\ and x = X3 
is increased, or at least not decreased. We know that U' does not intersect the 
spherical caps of the torus bubble, since they do not meet the region between 
the planes x = x\ and x = X3. Since the volume under U' is equal to the 
volume of the ball region of the torus bubble B, the volume in the remaining 
(disconnected) region must be larger than that in the original torus region T. 
There is no increase in area, contradicting Corollary 2.8. 



Figure 9. Reflection of a piece of T» could lead to extra intersections. 

Note that have used here the fact that two distinct Delaunay curves in- 
tersect locally at finitely many points. 

Assume now that the unduloid is the outer surface T . We again reflect 
it to get a new boundary surface U'. The total volume underneath U' is the 
same as that underneath U. Possibly Tj intersects the reflected surface U'. 
In this case the total volume underneath the union of the two surfaces has 
increased. Divide the region underneath U' so that the region under Tj is 
allocated entirely to B, as before, and all other regions are allocated to T. 
Then the total volume of B is preserved, that of T is increased, and the area 
is unchanged, a contradiction as before. □ 

We note that a somewhat similar argument was used by Athanassenas [7] 
in studying the stability of minimizing capillary surfaces. See also [41]. Our 
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situation has some additional complications due to the possible intersections 
resulting from the reflection. 

LEMMA 4.8. If a torus bubble has 1 < 60° and 6 2 < 60° then it is not 
minimizing. In particular, 62 > 60° in a minimizing torus bubble. 

Proof. If both 81 and 62 are less than 60° , the outer bubble T Q is generated 
by a nodary that becomes vertical twice. Proposition 4.6 implies it is unstable. 
Since 62 is the larger angle, 62 > 60°. □ 

Lemma 4.9. Suppose that T is a torus bubble with 9 2 > 120°. The the 
torus component is completely contained inside the smaller region of the double 
bubble whose larger outer sphere subtends an angle equal to 62- Similarly the 
ball component contains in its interior the larger region of this double bubble. 

Proof. Consider a double bubble which is symmetric around the x-axis 
and whose larger spherical cap S2 is on the right and subtends an angle of 
2 with the x-axis. Denote the smaller (leftmost) spherical cap by S Q and 
the spherical cap forming the interface by Si. Then S D is tangent to r G at 
(x\,yi). We will show that a torus region must be contained inside the smaller 
component B\ of this double bubble, the component bounded by S2 and Si. 
Otherwise, we will show that Tj and r G cannot rejoin after leaving (xi,yi), and 
thus cannot generate a torus region. See Figures 10 and 8. 

Given a vector in the upper half-plane, there is a unique circle perpendic- 
ular to the x-axis which is tangent to that vector. A Delaunay curve tangent to 
the same vector and generating a Delaunay surface with greater mean curvature 
than the corresponding sphere is a nodary, by Lemma 3.2. A Delaunay curve 
tangent to the same vector and generating a Delaunay surface with smaller, 
but still positive, mean curvature is an undulary or horizontal line. 

If the mean curvature of T is equal to that of S Q , then the mean curvature 
of Tj is equal to that of Si by Lemma 4.3. It follows from Lemma 3.2 that 
T Q , Ti and S D , Si coincide, and that there is no torus component. 

Suppose now that the mean curvature of T Q is smaller than that of S Q . 
Then the mean curvature of T is smaller than that of Si by Lemma 4.3. Since 
h Q > 0, also by Lemma 4.3, it follows that T Q cannot be a nodoid, and therefore 
t is a graph over the x-axis. In particular, r G does not recross the spherical 
cap S Q . The generating curve Tj for Tj leaves the point (xi,yi) on the positive 
(right) side of Si, and lies above S*j. Since S 2 forms part of the torus bubble, 
Tj must remeet r without crossing S 2 or itself if it is to generate a torus 
component with r G . If it is a graph over the x-axis then it cannot remeet r 
at all. If it is not a graph, then it is a nodary which generates a nodoid with 
negative mean curvature. Such a curve crosses S 2 before it can remeet r D . 
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Finally, suppose that the mean curvature of T Q is greater than that of S . 
Then Lemma 4.3 implies that the mean curvature of Tj is greater than that 
of Si, and Lemma 3.2 implies that both T Q and Tj are subsurfaces of nodoids 
and that near {x\,y\), T Q lies below S and Tj lies below Si. See Figure 11. 
The nodary r c generating T Q either crosses itself before leaving B\ or leaves B\ 
through Si. Similarly, the nodary n generating Tj either crosses itself before 
leaving B\ or leaves B\ through S Q . The curves r Q and Tj cannot self-intersect 
before remeeting, since they generate a torus component. Therefore either r G 
and Tj are contained in Si or they remeet for the first time after r G leaves B\ 
through S{ and Tj leaves B\ through S Q . The latter case is impossible, since 
Tj must cross r D before it can leave B\ through Si. It follows that the torus 
component is contained inside B\. 

The spherical caps Si in the double bubble and S\ in the torus bubble are 
disjoint, since S\ has boundary circle lying to the left of Si, and each is convex 
in opposite directions. So no part of the torus bubble intersects the interior 
of the larger, rightmost component of the corresponding double bubble, which 
must then be contained in the ball component of the torus bubble. This proves 
the last assertion of the lemma. □ 

Corollary 4.10. In a minimizing equal-volume torus bubble, 82 < 120°. 




Figure 10. If h is too small, a torus component cannot exist. 




Figure 11. The torus component is trapped inside a component of an 
associated double bubble if 82 > 120°. 

Proof. If 62 > 120° then Lemma 4.9 implies that the volume of the torus 
component is strictly less than the volume of the smaller component of the 
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double bubble whose spherical caps subtend angles 62 and (240 — #2)°- The 
other component of the torus bubble contains in its interior the larger com- 
ponent of this double bubble, which in turn contains more than half the total 
volume of the double bubble. It follows that the torus component has strictly 
less volume than the ball component of the torus bubble, contradicting the 
assumption of equal volume. □ 

Recall that 62 is the larger of the two angles subtended by the spherical 
caps, so the Corollary implies that both of 6\ and 62 are less than 120°. 

Lemma 4.11. In a minimizing equal-volume torus bubble, h Q > 0. 

Proof. We know from Lemma 4.3 that h > 0. If h = then r D is a 
catenary. Since ?/2 > V\ the catenary is given by a strictly increasing graph for 
{x > X2}, implying in particular that 62 > 120° which violates Lemma 4.9. □ 

Knowledge of h , 61,62 can be used to give a fairly accurate qualitative 
picture of the torus bubble. 

Proposition 4.12. A minimizing equal-volume torus bubble has the 
following properties: 

1. Tj has negative slope at (xi,yi) if and only if 61 > 30°. Its angle with the 
positive x-axis is 30° — 6\. 

2. n has positive slope at (0:2,2/2) if and only if 62 > 30°. Its angle with the 
positive x-axis is 62 — 30°. 

3. r has positive slope at (xi,yi) if and only if 61 > 60°. Its angle with the 
positive x-axis is 150° — 6\. 

4. t has angle at (2:2,2/2) with the positive x-axis equal to 62 — 150°. 

5. Tj is a graph. It has a unique local minimum if and only if 6\ > 30°. It 
never has a local maximum. If Ti is a nodoid then hi < 0. 

6. r always has a unique local maximum. It has a vertical tangent on the 
left if and only if 6\ < 60° . It never has a vertical tangent on the right. 
It never has a local minimum. 

7. If 6\ > 30° the local minimum y m \ n of Ti has value 

-h 

ymin ~ l + yr+W 

where fi is the force associated to Ti. 
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8. The local maximum 2/ max of t q has value 

_ l + yT + / ft 

2/max — , j 

K 

where f a is the force associated to T Q . 

Proof. The first four assertions follow from the fact that the surfaces form- 
ing a torus bubble meet at 120° angles along a triple curve. 

If n is not a graph then 9 2 > 120°, violating Corollary 4.10. If 6 1 > 30° 
then Tj has nonpositive slope at (x±,yi). Since 9 2 > Q\ the curve must have a 
local minimum in this case. 

Suppose Ti has a local maximum. Since 6 2 > 60°, Tj has positive slope 
at (x 2 ,y 2 ) and must have a local minimum as well. However 9 2 > 9\ so Tj 
must be longer than a period, contradicting Proposition 4.7. So it has no local 
maximum. 

If hi > and T is a nodoid, then we can compare Tj and T Q to a double 
bubble whose spherical interface Sj coincides with S±, and apply Lemma 3.2 
as in Lemma 4.9. T is trapped inside the component B 2 to the right of the 
interface and T Q is trapped inside the component B\ to the left of the interface, 
as in Figure 12. In particular, h Q < in this contradiction. So if T is 

a nodoid then hi < 0. 




Figure 12. A double bubble traps T and T D if T is a nodoid with hi > 0. 

If 6\ < 30° then Tj has positive slope at (x\, y±). If it has a minimum then 
it must pass through a maximum first, so it has no minimum. 

r is always oriented upwards at and downwards at (x 2 ,y 2 ), so it 

always has a local maximum. If it has a local minimum then it would contain 
two local maxima, violating Proposition 4.7. If 9\ < 60° then the curve starts 
out to the left, therefore it has a vertical tangent at which it changes direction 
from left to right. Proposition 4.6 implies that there cannot be a second vertical 
tangent, so it has no vertical tangent on the right, as claimed. If 9\ > 60° then 
r initially is going right. If it has a vertical tangent on the right then the 
angle 9 2 must be less than 60°, contradicting the assumption that 9\ < 9 2 . 
This proves the assertion concerning vertical tangents of t q . 
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To calculate the local minimum y m j n of , we use equation 2 of Section 3 
for the force of Tj, 

fi = hiy 2 - 2ycosa . 
At a minimum, a = and cos a = 1, so 

/i = ^i2/ 2 - 2y. 
If /i 7^ this quadratic expression gives two roots for y: 



1 ± y/TThiTi 

v = — i • 

The value of the minimum is given by 

(12) WTTM_ 

Hi 

whatever the sign of hi. If hifi > 0, as Corollary 3.3 implies happens at a local 
minimum for a nodoid, then the other root is negative and meaningless. If 
hifi < 0, as Corollary 3.3 implies happens at a local minimum for an unduloid, 
then the other root corresponds to the local maximum of the unduloid. A 
double root occurs only in the case of a cylinder, equation 12 becomes linear 
when hi = 0. An equivalent expression for the minimum, which holds also 
when hi = 0, is 

(13) ymhl = l + VTTWi' 

For a local maximum, a similar analysis gives that y max is one of the two roots 

(14) 9 _ = i±4±M. 

Since by Corollary 4.19 T is always a nodoid, h a > and f a > 0, there is 
a unique positive root, and 



_ l + yT+ f h 

2/max — , • I— I 

h Q 

Proposition 4.13. Given 6\ and h a , there is at most one corresponding 
minimizing equal-volume torus bubble. 



Proof. Given 9\ and h Q , hi = 2 — hi is also determined. Thus S\,T and 
Tj are uniquely determined. The spherical cap S2 is uniquely determined by 
the second intersection point (x2,U2) of r D and Tj, assuming that this second 
point exists. If r D and Tj do not meet in a second point then there is no torus 
bubble corresponding to 9\ and h Q . In general, r D and Tj may intersect in many 
points, and it is necessary to show that only one of these can possibly lead to 
a torus bubble. 
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If t and Tj are graphs, then it is clear that their first intersection point 
with x coordinate larger than x\ is the unique intersection point which can give 
a torus bubble. In general, Tj is a graph, but not t . However Proposition 4.12 
shows that r has no vertical tangents for x > xi, so that both curves are 
graphs in this region and the same argument applies. □ 



Figure 13. Some eliminated possibilities for cross-sections of torus bubbles. 

Proposition 4.14. // a torus bubble has 6i+6 2 > 180°, then it is not a 
minimizer. 

Proof. Under this assumption there is a value of y realized by three dis- 
tinct points on the circular arcs generating the spherical caps. Since 62 > 6\ we 
have in particular that 62 > 90°. We know from Lemma 4.9 that 62 < 120°. So 
90° < 6» 2 < 120°, 60° < 0i < 120°. Then Lemma 4.12 implies that the curves 
Ti and t are both graphs. If 61 + 62 > 180° then there are three points on 
S\ U S 2 with y-value equal to y 2 , (ai,y 2 ), (02,2/2), (03,2/2), with ai < a 2 < a 3 
and d2 = x 2 . Reflect the part of the bubble between the planes x = a\ and 
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x = ci2 through the plane x = (x±+X2)/2 to get a new piecewise-smooth surface 
with the same area enclosing the same volumes, as in the proof of Lemma 4.7, 
but without needing to worry here about introducing self-intersections. The 
new piecewise-smooth surface has three surfaces meeting at an angle not equal 
to 120° and so is not minimizing, a contradiction. □ 

Propositions 4.6, 4.7, 4.14, 4.12, Lemma 4.8 and Corollary 4.10 combine 
to rule out many possibilities for the shapes of minimizing torus bubbles. See 
Figure 13 for cross-sections of some of these eliminated possibilities. 

Lemma 4.15. For any constant 6, the function /(c) = 6(1 — (?) + 
c^/3(l — c 2 ) is at most two-to-one on (—1, 1). 

Proof, /'(c) = —26c — \/3g(c), where 

g(c) = (2c 2 - 1)/VT^. 

Since g"(c) = 3(1 — c 2 ) -5 / 2 > 0, g is convex on (—1, 1). 

A line intersects a convex curve at most twice, so /' has at most two roots 
on (—1,1). Thus / has at most one local maximum and one local minimum 
on (—1, 1). Combined with /(0) = /(l) = 0, this completes the proof. □ 

We calculate in the next proposition the value of the forces f a of T Q and 
Si of Ti. 

Proposition 4.16. For the torus bubble determined by 0\ and h Q , the 
forces f of T and fi ofTi are given by: 

So = (K - 1)(1 - Cl 2 ) + c lv /3(l- Cl 2 ) 

and 

fi = {hi - 1)(1 - ci 2 ) - 01^3(1 -ci 2 ) = -So 
where c\ = cos 6\ . 

Proof. We first calculate the force f a of T a . At (xi,yi), y\ = sin#i and 
a = 5"7r/6 — 0\ is the angle between the positive x-axis and the graph of r D , 
therefore 

So = Kyi 2 - 2yicosa 

= h D sin 2 9\ — 2 sin#i cos(57r/6 — 6\) 

= h a (l — cos 2 6\) - 2 sin 6*i ({-V3/2) cos 0i - (1/2) sin (-^i)) 

= h (l - ci 2 ) + ci^3(l - ci 2 ) - (1 - ci 2 ) 

= (/ io -l)(l-ci 2 )+ci V /3(l-ci 2 ) . 
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Similarly, 

fi = hiyi 2 - 2y 1 cos a 

= (^-l)(l- Cl 2 )- Cl ^3(l- Cl 2) . 

Since h a — 1 = 1 — hi = — (hi — 1) by Lemma 4.3, f = —fi. □ 

Note. The fact that f = —fi is an example of a general "balancing 
principle" for constant mean curvature surfaces; see [25]. 

Lemma 4.17. For the torus bubble determined by 6\ and h , the angle 2 
is one of the (at most) two solutions of the equation 

/(cos 01 ) = /(cos0 2 ), 

where 

/(c) = (h Q - 1)(1 - c 2 ) + Cv /3(l-c 2 ). 



Proof. We calculate the force f a of T Q at both (x±,yi) and (x 2 ,y 2 ). 
At (xi,yi), 

/ = - 1)(1 - Cl 2 ) + Cl ^3(l- Cl 2 ), 
where j/i = sin 0i and C\ = cos 0i . A similar calculation at 62 shows 

So = (K - 1)(1 - c 2 2 ) + c lv /3(l-c 2 2 ). 

Setting the two forces equal and applying Lemma 4.15 gives that there are at 
most two possible values for c 2 = cos0 2 - Since < 2 < 180°, where the cosine 
function is monotonically decreasing, this implies that there are at most two 
possible values for 2 . □ 

One solution occurs when the angles are equal. An example of such a 
symmetric torus bubble, depicted in Figure 7, occurs with 9\ = 62 = 90°, and 
h a approximately 1.9848. However these symmetric solutions never lead to a 
minimizing torus bubble by Corollary 4.5. 

Recall that we remarked after Corollary 4.5 that we can without loss of 
generality restrict attention to nonsymmetric torus bubbles. This means that 
we are working in the case where 0i < 2 . For each 0i, h Q , there is at most 
one value of 62 which can potentially lead to a stable torus bubble. 

PROPOSITION 4.18. // a minimizing torus bubble has 60° < 0i < 120° 
then h Q > 1 — \/3cot 0i and hi < 1 + y/3 cot 0i. The same statement holds with 
62 replacing 9±. 



490 



JOEL HASS AND ROGER SCHLAFLY 



Proof. We again apply Lemma 3.2 to compare the torus bubble with a 
double bubble whose smaller spherical cap S a subtends an angle of 0\, as in 
the proof 1 Lemma 4.9. See Figure 14. If hi is larger than the mean curvature 
of Si, then Tj is a nodoid contained inside the smaller component of the double 
bubble, as shown in Lemma 3.2. See Figure 14. It follows that r G cannot 
remeet t\ and the two curves cannot generate a torus component. □ 




Figure 14. The mean curvature of T Q cannot be too small if r Q and Tj 
are to meet again. 

Corollary 4.19. T Q is a nodoid. 

Proof. If 60° < #i < 120° then Proposition 4.18 shows that the mean 
curvature of T Q is greater than that of the sphere which is tangent to it at 
(x±,yi), so Lemma 3.2 implies that T Q is a nodoid. If Q\ < 60° then r D starts 
out heading to the left. It must become vertical, hence T Q must be a nodoid. 

□ 

Lemma 4.20. Let N be a subcurve of a nodary generating a nodoid of 
mean curvature h > 0, running between a vertical tangency on the left at 
(xv,Vv) and o, vertical tangency on the right at (x' v ,y' v ), with a unique local 
maximum at (xm,Um)- Then a circle C of curvature h which is tangent to 
N at (x v ,y v ) lies underneath N. Moreover x v — x m is greater than the radius 
1/h of the circle. 

Proof. In a deleted neighborhood of (x v ,y v ), the curvature of N is less 
than that of C, so that C locally lies beneath N. If C crosses N between x v 
and x' v , let C be the subcurve of C starting at (x v ,y v ) and running to the first 
point of intersection of C with N. Let C[ = C + (0, t) be a vertical translate 
of C by a distance t along the y-axis. Let T = sup{t : C[ n N / 0}. Thus C' T 
is tangent to N at a point P and lies above N near P. The curvature of N 
at P is given by k m = h — k p < h. Since the curvature of C' T is greater than 
that of N, it cannot lie above it, and it follows that X ]\/f Xy IS larger than the 
radius of C, i.e., xm — x v > 1/h. □ 
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The next proposition rules out the existence of minimizing torus bubbles 
with h Q and Q\ both close to 0. The numerical calculations that we will apply 
are badly behaved in this region, so we use this geometric argument to exclude 
it. 

PROPOSITION 4.21. // a torus bubble has h < 0.2 and 9\ < 5.7°, then it 
is not a minimizer. 



Proof. We will show that before Tj and r G rejoin at (x2, 3/2)1 either r c has 
two vertical tangencies, or r, traverses more than a full period. Note that our 
assumption implies that hi > 1.8 

Proposition 3.4 implies that T{ is a graph and its period is no longer than 
lit I hi < 3.5. The bubble is not a minimizer if X2 — x\ > 3.5 by Proposition 4.7. 
The initial angle of Tj from the positive x-axis is 30° — 0i, so it cannot turn 
(anti-clockwise) through an angle of more than 6O°+0i. When they meet again 
Ti and t intersect at an angle of 120°. It follows that r does not intersect Tj 
before it turns (clockwise) through an angle greater than 180° — 0±. Therefore 
the two curves cannot intersect before the nodary r G reaches its maximum, 
which happens when r has turned through an angle of 150° — 9±. Let (xm,Vm) 
be the point where r reaches its maximum and (x v ,y v ) the point where r first 
goes vertical. Then xi — x\ > xm — x\. 

Now h a = k p + k m where k m is the curvature of the generating curve in 
the xy plane and k p is the principal curvature due to rotation around the axis. 
At (x v ,y v ), k p = and h Q = k m < .2 by hypothesis. 

Lemma 4.20 implies that xm — x v > l/h Q > 5. 

The assumptions h a < .2 and 6\ < 5.7° give an upper bound for the force 
fo of T Q . Calculating at height y\ = sin^i < 6\ < 5.77r/180 < .1 gives 



\fo\ = |yi 2 h o -2/icos(3O o + 0i)| 

< \yi 2 K - yiV3/2| 

< {.l) 2 h Q + .1V3/2 

< .09. 

We bound the value of y v from above by solving f Q = y^h . 

y 2 v = f /ho < -09 /ho, 

so 

y v < y/.09/h o . 

We next bound from above by noting that t q is convex between x v 

and x\ by Lemma 3.4. This implies that (x v ,y v ) lies to the right of the point 



492 



JOEL HASS AND ROGER SCHLAFLY 



(x,y v ) at which the line of slope f (x\) has height y v , so that 

\x v -xi\<\x-xi\ = \(y v - y\)/f {xi)\ < V3(y v - yi) 
< V3y v < V3^/.09/h o = yf.27/h . 

This implies that 

xm ~ xi = (x M ~ x v ) - (xi - x v ) > l/h - ^.27/h Q . 

For < h Q < .2, differentiation shows that the function l/h Q — \J .27 jh Q 
is decreasing and its minimum value occurs at h Q = .2, so that xm — x\ > 
1/(0.2) - ^.27/(0.2) > 3.8. Thus r travels a distance of at least 3.8 to the 
right before hitting Tj. But the period of Tj is less than 3.5. We have shown 
that either r D changes direction and has two vertical tangencies or Tj contains 
a full period. In the first case it is not a minimizer by Lemma 4.8, and in the 
second case it is not a minimizer by Proposition 4.7. □ 

The next lemma will be used to give bounds on the size of a circle enclosing 
a loop of a nodary. 

LEMMA 4.22. // an embedded curve 7 in the positive quadrant of the 
plane is tangent to the y-axis at the origin and has curvature k > 1/r > 0, 
then it does not meet C(r), the circle of radius r centered at (0, r), at any point 
other than the origin. 

Proof. If it does cross C(r), denote by 7' the subcurve of 7 starting 
at the origin and running to the first point of intersection with C(r), as in 
Figure 15. The subcurve 7' is a graph over the x-axis since it is embedded. 
Let 7t = 7' + (0, t) be a vertical translate of 7' by a vertical distance of t. 
Let T = sup{t : i t n C(r) / 0}. Then j' T meets C(r) at a point P in its 
interior, but not the interior of C(r), and thus is tangent to C(r) at P. Since 
the curvature of j' T is greater than that of C(r), this is a contradiction. □ 




Figure 15. The curvature of the curve 7 cannot be greater than 1/r 
if it crosses C(r) twice. 

Lemma 4.23. // r G has curvature greater than k, and T{ has curvature 
less than —k, then their union is contained in the intersection of two circles of 
radius r = 1/k meeting at 120°. 
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Proof. Construct circles C\ and C2 of radius l/k, meeting at 120° at 
(xi,yi), with C\ tangent to r G and C2 tangent to Tj as in Figure 16. Let D 
denote the line segment joining the two intersection points of C\ and C2. By 
Lemma 4.22 we know that each of r G and Tj cannot leave the intersection of 
these two circles without first crossing D. Since each is convex in C\ n C 2 it 
follows that they remeet before either one leaves C\ f~l C 2 . It follows that the 
intersection of the torus component with the upper half-plane is contained in 

Cinc 2 . □ 




Figure 16. The torus component generating curves and r are 
trapped inside two circles meeting at 120°. 

Lemma 4.24. If a torus bubble has hi < — k < and the intersection of 
the torus component with the upper half -plane has minimum y-value equal to 
yo with 2/0 > l/k, then the torus component contains volume 

v <2.5ir(y + —- -)(- 



2 ky - 1 ky - 1 ' 



Proof. In a torus bubble with hi < — k < we have that \h a \ = —hi + 
2 > k. For a Delaunay curve passing through a point (x, y) in the upper 
half-plane, 

\k m \ ^ \h\ — \k p \ 
= \h\ 



y^l + ij 2 

> 

For the curves T{ and r G we have y > yo, so that 
(15) \ km \>k--. 

yo 

So the curves r G and Tj each have curvature with absolute value greater than 
(ky - l)/y . 
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By Lemma 4.23 the intersection of the torus component with the up- 
per half-plane is contained in the intersection w of two circles of radius r = 
Uo/{kyo — 1) meeting at an angle of 120°. The area of w is 



= (2tt/3 - V3/2)r 2 < 1.25(- 



yo 



The diameter of w is 

kyo - 1 

Since w is centrally symmetric and has minimal y- value smaller than or equal 
to yo, the center of mass of w has y-coordinate at most yo + d/2. Pappus' 
theorem for the volume of a solid of revolution obtained by rotating this area 
around the x-axis gives that 

• < 2, feo + |)1.25(^ T f, 

< 2.5,( W + ^ »)(»»)» 
2 kyo-1 ky - 1 

We apply these volume estimates to reduce the range of possible mean 
curvatures in a minimizing equal-volume torus bubble. 

Proposition 4.25. In a minimizing equal-volume torus bubble, h Q < 10. 



Proof. The hypothesis h Q < 10 is equivalent to hi > —8. Suppose to the 
contrary that hi < — 8 in a minimizing equal- volume torus bubble. 

Proposition 4.12 states that Tj is a graph, so that the volume Vi of the 
ball component B is strictly larger than the volume of the region under S*2. 
Lemmas 4.9 and 4.8 imply that 60° < 2 < 120°. It follows that v { > 5vr/24, 
the volume under a unit radius spherical cap subtending an angle of 60°. We 
now estimate the volume v of the torus component. The curves Tj and r G 
generating the torus component contain the point (x2,y 2 ) where y/3/2 < y 2 < 
1. Their maximum and minimum y-values differ by less than 2/\hi\ < 0.25 by 
Proposition 3.4, and therefore both curves are contained in {y > y/3/2 — 0.25}. 

The expression - — — — is decreasing with respect to yo when yo > 1/k, and 
kyo-1 

decreasing with respect to k, and also 0.6 < y/3/2— 0.25. We apply Lemma 4.24 
with k > 8 and yo > 0.6, giving that 

v < 2.57r(y + — -, r)(- 



2 kyo-1 kyo-1' 

n , „ v 7 ^ 0.6 w 0.6 x2 
< 2.5vr(l + ^— — -)(- 



2 4.8 - rM.8 - 1' 
< 0.08vr . 
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Thus v < 0.087T < 57r/24 < Vi, contradicting the assumption that the 
torus bubble encloses equal volumes. □ 

Our algorithm will need the following rather technical result about the 
slope of a generating curve. 

Lemma 4.26. Suppose that r% satisfies hi < 2 and < 9\ < 4°. Let 
a be the angle made by Tj with the positive x-axis and let a.\ be the value of 
this angle at (x\,yi). Then a > a.\ for as long as Tj remains below height 
{y = 1/2}. In particular, a > 26° below this height. 

Proof. First note that if hi < then Tj is a nodary or catenary, and the 
angle made by with the positive x-axis is increasing as one moves along Tj, 
which implies the conclusion of the lemma. So we can assume that < hi < 2. 

Let y be any point on Tj. Then 

hiy 2 — 2y cos a — fi = . 

A given value of cos a is realized by at most two y values solving this quadratic 
equation. Suppose that as the T{ goes from a minimum to a maximum, the 
angle 26° < ct\ = 30° — 6\ < 30° occurs twice. The two y-values which share 
the same ct\ both solve the above equation. Thus they sum to 2cosai//ij > 
2(.8)/2 > .8. Since the smaller root has y-coordinate sin(4°) < .1 the larger 
root is greater than .7 > 1/2. The angle a is increasing near the smaller root. 
Thus Ti subtends an angle of at least ct\ with the x-axis as long as T{ remains 
below height 1/2. Since ct\ = 30° — 9\> 26°, the final conclusion holds. □ 

5. Computation 

We have shown that to solve Conjecture 3 in the equal volume case it suf- 
fices to show that certain torus bubbles are less efficient at enclosing two equal 
volumes then a double bubble. In this section we will show how these torus 
bubbles are ruled out by a computation. The algorithm for this computation 
is given below. 

The basic idea is to consider a domain of torus bubbles corresponding to 
a product of small intervals in each of 9\ and h Q , and to calculate as much 
as possible about the geometry of the corresponding torus bubbles. Various 
calculations are then applied which rule these torus bubbles out as potential 
minimizers. The accuracy of these calculations depends on the size of the 
domain rectangle we start with. The computational scheme will succeed if 
these can be chosen small enough to get sufficient accuracy, yet large enough 
that a reasonable number of them cover all the possibilities. 

The computation was performed using double precision floating point 
numbers. The fundamental data type used is the IEEE 754 32-bit real number, 
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see [5] . This is a binary representation with a 23-bit mantissa (plus an implied 
leading bit), an 8-bit binary exponent, and a sign bit. The (finitely many) 
real numbers that can be represented in this standard with no error are called 
representable. The IEEE standard specifies that the add, subtract, multiply, 
divide, and square root operations be performed as if done exactly and rounded 
to the nearby representable number according to the rounding mode in effect. 
There are three rounding modes that can be chosen: up, down, or nearest. The 
values +00 and —00 are representable and behave in a specified way. Floating 
point exceptions are masked, but flags are sticky and available for clearing and 
inspection. Most computers in use today implement the IEEE standard. We 
denote the result of rounding a real number x down to a representable number 
smaller or equal to x by x, and similarly denote by x the result of rounding x 
up to a representable number greater or equal to x. 

Combined with the methods of interval arithmetic, see Moore [29] and 
Alefeld-Herzberger [2], the IEEE standard allows numerical calculation with 
exact bounds on accuracy. Interval arithmetic is a method by which a real- 
valued function on the reals can be extended to an interval-valued function of 
intervals. An interval is formed from two representable reals. Mathematically, 
it represents the closed interval between the two reals. The add, subtract, mul- 
tiply, divide, and square root functions are extended to intervals by the IEEE 
operations on reals along with directed rounding. We only use the round-to- 
nearest mode in calculating averages, (see the procedure avgwt below), as these 
are only used to divide intervals and do not need to have bounds calculated. 

Generally, we say that an interval X is an extension of a real number x 
if x € X and an interval-valued function of intervals F is an extension of / if 
xGl4 f(x) G F(X). 

In this section, real numbers are denoted by lower case, and intervals with 
representable real endpoints by upper case. The reals embed into the set of 
intervals by mapping a real to the smallest interval containing it. The lower 
and upper bounds to an interval are denoted with lower and upper bars, so 
X = [X_, X]. Arithmetic operations on representable reals and intervals are in- 
terpreted according to IEEE and interval rules, not by the usual mathematical 
definitions. 

As an example, the sum of two intervals A = [A, A] and B = [B_, B] is 
given by 

A + B=[ A + B ,A + B], 

where the left endpoints are rounded down when added, and the right endpoints 
are rounded up when added. 

Other operations on reals are extended to intervals as in [29] and [2]. 
Relations are interpreted positively, so for example when X and Y are intervals, 
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X < Y means that for any x 6 X and any y £ Y , x < y, and 1/7 means 
that for any x £ X and any y £ Y, x / y. 

We have equivalent expressions X < Y 44> X < Y_, and I / 7 ^ 
IflF = 0. The union of intervals X and Y is somewhat nonstandard; it is 
the smallest interval containing both sets, and is denoted XUY. The interval- 
valued functions Absolute Value(X), Max(J, 1"), M.in(X, Y), and the 
intersection of two intervals X (~)Y are defined in the standard way, without 
rounding. There is no IEEE standard for transcendental functions, so we 
designed our program to avoid all calls to trigonometric functions. 

The operations of the IEEE floating point standard can lead to undefined 
operations. For example the quotient 0/0 and the product x +oo result in 
an output of NAN (Not A Number), and are signaled by the presence of an 
exception flag. Operations of interval arithmetic can also lead to questionable 
operations, for example when the quotient of two intervals A/B is calculated 
and B contains 0. Our implementation of the division operator on intervals 
returns [— cx),+oo] if B contains 0. The algorithm draws no conclusions when 
such a division occurs, but rather calls for a subdivision of the input into inter- 
vals of smaller size, where the operations are repeated with greater accuracy. 
Exceptions such as overflow and underflow, sometimes an issue in computer 
assisted proofs, are not an issue in our algorithm. In any case they do not 
occur. 

Bounds obtained from repeated applications of interval arithmetic are 
potentially far from sharp, especially when using wide input intervals. If just 
a single value needs to be computed rigorously, one can work with intervals 
whose width is comparable to the roundoff error of a single operation. For 
our purposes these types of intervals are far too thin, and would make the 
number of calculations we need impractical. Our intervals are "fat" in the 
terminology of Fefferman [15], meaning that they are sized by the scope of the 
problem we are solving rather than the size of the computational rounding. 
Often the intervals are wide enough to make a perfectly good formula look like 
nonsense. Sometimes it is possible to narrow the intervals under consideration 
because we have knowledge about what the possible legitimate values are that 
could arise during a computation. For example, consider a formula involving 
square roots. Assuming we're not using imaginary numbers, the validity of 
the formula presupposes that the argument to the square root is nonnegative. 
However, when we pass to an interval extension of the formula, the interval 
argument to the square root will often include negative values, even though 
these values cannot arise from the problem we are interested in. The square 
root of a negative number would normally trigger an exception, but the nature 
of the calculation allows us to define the interval square root function Sqrt to 
discard any negative portion of an interval argument. The justification for this 
apparent sleight-of-hand rests on the validity of the original formula over the 
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reals. If our theorems tell us that the quantity whose square root we are taking 
is nonnegative, we are justified in truncating the interval to exclude negatives. 

The interval function which always returns [— oo,+oo] can be used to 
represent any function, but not very usefully. We use this interval to return 
the value of a division by an interval containing zero. An interval [a, b] is empty 
if a > b. All of our functions take nonempty interval arguments and return 
nonempty interval results except for the intersection operation, which returns 
an interval which can be empty. 

The value of the nonrepresentable constant y/3 in our algorithm is ex- 
panded to a narrow interval. Other constants used are represent able. Thus we 
use the expression 5H a < 1 rather than H Q < 0.2, which involves the nonrepre- 
sentable constant 0.2. Representable constants such as 5 behave identically in 
interval arithmetic to intervals where both endpoints have the same value, such 
as [5,5], and operations involving both intervals and representable constants 
are treated as if the constant was first converted to a one-point interval. 

The verify statement is used to assure that a particular condition holds. 
If the condition fails, then the entire program is stopped (aborted). If no such 
failure occurs, the presence of such a verify statement constitutes a proof that 
the condition holds for the given inputs. One use of the verify statement is 
in the procedure avgwt defined below, which gives a weighted average of two 
intervals. To apply this we want to be sure that, as we expect in the cases 
where it is applied, the first interval is strictly smaller than the second. We 
check that this is so with a verify statement. Running our program results in 
no violations of any statement which is tested with a verify statement. 

The program examines the set of all torus bubbles to see if any can be a 
minimizer. A range of hypothetical torus bubbles is specified by intervals ©i 
and H , according to Proposition 4.13. To avoid unnecessary use of trigono- 
metric functions our algorithm uses an equivalent parameterization of the space 
of torus bubbles by intervals Y\ and H Q , where Y\ = sin@i. Rather than solve 
the mean curvature differential equation directly to find intervals (X2, Y2) con- 
taining the point (£2,2/2) where Tj and r G intersect, we assume existence of a 
torus bubble, derive Y2, and then deduce X2 from numerical integrals. Volume 
calculations for Vi and V are obtained from additional numerical integrals. 
At each stage in the calculation, we check whether the torus bubbles can be 
rejected based on the instability results of Section 4. If not, we test whether 
the volumes are equal. As it turns out all torus bubbles are rejected for these 
reasons, so area calculations are not necessary and are not present in the algo- 
rithm we present. Most torus bubbles are rejected either because ©2 is out of 
range, or because the x-displacement of Tj and r D differ, or because the volumes 
enclosed in the two regions differ. 

The computation of ©2 and Y2 is based on an analysis of the force function. 
The integration process used to calculate the value of X2 is complicated by the 



DOUBLE BUBBLES MINIMIZE 



499 



fact that the curves may not be graphs, so we have to choose an appropriate 
parameterization. The most convenient parameterization is in terms of y, 
because the ODEs satisfied by Delaunay curves involve only y, and because 
the y-coordinates of the endpoints are specified by Gi and O2. This allows 
x-displacement and volume to be expressed directly as integrals in terms of y, 
making it unnecessary to actually generate an ODE solution which produces 
Tj and t . 

The curves Tj and r are generally not graphs as functions of y, so the 
integrals we use have singularities at points where y = 0. A change of variables 
resolves the singularity, as worked out in Proposition 3.7. 

Another difficulty is that the geometry sometimes degenerates at the 
boundary of the regions we are examining. For example, we need to exclude 
torus bubbles with 6\ arbitrarily close to 0, but some of our formulas become 
singular when 0\ = 0, and don't make sense there. If the interval 0i contains 
0, then various intervals representing y-coordinates along the Delaunay curves 
will also contain 0, and the integrals give results which will not allow torus 
bubbles to be eliminated. We get around this problem by using crude esti- 
mates which apply in a small zone near the x-axis. Near the x-axis we have 
some information about the initial angle and the concavity of the Delaunay 
curves, and this suffices to give crude bounds to their behavior near the axis. 
We compute integrals only when at an appropriate distance above the x-axis. 
This approach suffices for all cases except for when h Q is also near 0. In that 
case we rely on Proposition 4.21. 

We now describe in pseudo-code the main interval procedures used in our 
algorithm. As noted above, aside from the procedures described in this section 
our algorithm makes use only of standard arithmetic operations on intervals, 
namely addition, subtraction, multiplication, division, absolute value, maxi- 
mum, minimum and square. 

The procedure avgwt calculates a number in between two inputs, used 
later to subdivide intervals. It returns a representable real number which is 
a weighted average. For convenience, we set the rounding mode to round to 
nearest in this procedure. In some applications of this procedure, an interval 
is assigned the output of the avgwt procedure, eg in Step 7 of the procedure 
DivideAndCheckRectangle. The resulting interval then consists of a single 
point. 

real procedure avgwt 
input : X , Y , w 

verify X < Y 
z:= (l-w)X + wY 
return z 



500 



JOEL HASS AND ROGER SCHLAFLY 



The interval procedure Sqrt(X) returns an interval containing {^fx : x £ X 
and x > 0}. The real procedure width calculates an upper bound for an 
interval's width. 

The interval procedure Compare checks whether the interval X is to the 
left of the interval Y, returning A if yes, B if X is to the right of Y, and the 
union AUB if the intervals overlap. 

interval procedure Compare 
input : X, Y, A, B 

if X < Y then return A 
else if X > Y then return B 
else return AUB 

Lemma 5.1. For any x,y,a,b, in any X,Y,A,B, respectively, if x < y 
then a € Compare(X, Y, A, B), and if x > y then b € Compare(X, Y, A, B). 

Proof. The lemma follows from a straightforward case by case analysis. □ 

Integrate is the basic method for numerical integrals. It gives relatively 
wide intervals, but is good enough for our purposes. Its input is an interval- 
valued function of one interval variable, and two intervals which give the limits 
of integration. 

interval procedure Integrate 
input : F,A,B 



return 0UF(A)Width(A) + 0UF(5)Width(5) + RH 

Proposition 5.2. For any integrable f and real numbers a, b, with a < b, 
and interval extensions F, A, B, the integral of f from a to b is contained in 
Integrate (F, A, B). 



R:=0 

H := (B-A)/32 
if (H > 0) then 



begin 

for i := to 31 do 



R:=R + F(A + (OUH) + iH) 



end 



Proof. If A < B, 




■b 



■A 



rB_ rb 

f dx + _ f dx + I f dx . 
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rA rb 

The term / / dx is contained in 0UF(v4)Width(^4) and similarly I f dx 

Ja Jb_ 
is contained in OUF (B) Width (B). Note that we need to replace F(A) and 

F(B) by 0U-F(^4) and OUF(B) to deal with the issue that a can lie anywhere 

in A and similarly b in B. Taking the union of F with ensures that our sum 

contains the intervals F([a, A]) and F([B_, b]). 

f- 

The value of /_ / dx is bounded by lower and upper Riemann sums, using 
Ja _ 
32 equally sized intervals. The term F{A + (0L)H) + iH) evaluates F on an 

interval which contains [A + ih, A + (i + l)h] for any h € H, and so contains 

both an upper and a lower bound for F on this interval. The RH term in the 

last line, when added to the previous terms, then gives an interval containing 

rb 

both lower and upper Riemann sums for / / dx. 

J a 

Otherwise, A < B is false, and since a < b, A and B must overlap, so we 
can decompose 

rb rc rb 

/ f dx = f dx + f dx 

J a J a J c 

where a < c < b and c G A n B. Since B_ < A, R = and the RH term in 
the last line vanishes. The result follows from [a,c] C [a, A], [c,b] C [JB, b] and 
inclusion monotonicity. □ 

The next several procedures define functions to be integrated. 

The procedure Dx is a straightforward interval extension of the expression 
for dx/dy given in equation 6 in Proposition 3.6. It will be used to find the x- 
displacement from {x\, y±) to (x2, V2) by applying equation 6 to both Ti and r G . 

Volumes under Delaunay surfaces are computed using the volume integrand 
Dv, which is obtained from Dx as in Proposition 3.6. It is convenient to have 
Dv calculate volume divided by it. Since volumes are only compared to one an- 
other, we can avoid the unnecessary step of multiplying by the nonrepresentable 
constant ir. The volume element computation required that dx/dy > for the 
volume to have a positive sign. The curve r G will sometimes start out in the 
negative x direction, when 6\ < 60°, and after passing through a vertical tan- 
gency switch to going in the positive x direction. Calculating the integral of 
Dv in this case will give the volume under the part of the nodoid generated 
by the arc going right minus the volume under the arc going left. This is just 
what we want to compute the volume V in the torus component in the case 
when such an overhang exists. 

At local minima and maxima, Dx and Dv are singular. Y m [ n and Y max are 
global variables which are accessed by Dx m i n , Dx max , Dv m i n and Dv max . So 
near local extrema, we use a change of variables to reparameterize the curve, 
as in Proposition 3.7. Procedure Dx m j n calculates dx/dz near a minimum 
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and procedure Dx max near a maximum, where z is given in Proposition 3.7. 
We do an algebraic manipulation to make the denominator in the square root 
more computationally effective, namely we replace the term HY — 2 + i7Y max in 
procedure Dx max by the equivalent expression H(Y + Outcrmin) • Here Yutermin 
represents the minimum y-value of the nodoid which has a maximum at Ym ax . 
The equivalence of the two expressions is derived using Proposition 3.4 as 
follows: 

HY-2 + HY max = -1 (Y + y max ) -2 = H(Y + Y outermin ). 

'max ^outcrmin 

The volume calculation procedures Dv mm , Dv max give the formula for the 
volume integrand used near a critical point. The volume calculations are ap- 
plied only to graphs over the x-axis. 

interval procedure Dx 
input : Y,H,F 

T := HY 2 - F 

return T/Sqrt((2Y + T)(2Y - T)) 



interval procedure Dx m j n 
input : Z,H,F 

Y '.= Y m i n + Z 
T := HY 2 - F 

return 2T/Sqrt((2Y + T)(2 - HY min - YH)) 



interval procedure Dx max 
input : Z,H,F 

Y . — Y max Z 
T := HY 2 - F 

return 2T/Sqrt((2Y + T)H(Y + F uter m in)) 



interval procedure Dv 
input : Y,H,F 

return Y 2 Vx(Y,H,F) 



interval procedure Dv m i n 
input : Z,H,F 

return (Y min + Z 2 ) 2 Dx min (Z,H,F) 
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interval procedure Dv max 
input : Z,H,F 

return (Y max - Z 2 ) 2 Dx max (Z, H, F) 

The next procedure is used for testing interval rectangles of torus bubbles to 
see if they are potential minimizers and rejecting these intervals if they can be 
shown to have a property which rules out any torus bubble contained in them. 
The input is a pair of intervals C\ and H Q , which give a range of possible values 
for cos 9\ and h Q , respectively. The interval of cosine values C\ is used, rather 
than an interval of Q\ values, to eliminate calls to trigonometric functions. A 
variety of tests are performed based on the analysis of previous sections. The 
various intervals do not necessarily give sharp estimates, so it is not a priori 
clear how many rectangles will be rejected. As it turns out, the tests are sharp 
enough to reject all rectangles. 

Boolean procedure Che ckRect angle 
input : C\ , H Q 

1 . if lOOOCi > 996 and 5H Q < 1 
then return REJECT 



2 . Hi : — 2 — H Q 

Yl := Sqrt(l - C 2 ) 

Ft := (Hi - l)Y 2 - CiYiSqrt(3) 

F := -Ft 

C 2 ■= [-dud] n [-.50.5] 

T := (2(Hi - m + 3)/(3 + (H t - l) 2 ) - (1 - d 2 ) 
if (C 2 2 + T) / 1 return REJECT 

3. Y 2 := Sqrt(Tn [0,1]) 




if 



C2 is empty then return REJECT 



4. 



if 



Ci < .5 and H a < l-Sqrt(3)Ci/Yi 
then return REJECT 



5. 



if 



Width(C 2 ) > .5 then return NDRESULT 



6. 



W ends := (1 - d) 2 (2 + d)/3 + (1 - C 2 ) 2 (2 + C 2 )/3 
Y min := -Fi/(1+ Sqrt(l + F l H l )) 
Y max := (1+Sqrt(l + F H ))/H 
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Y : = Compare(Ci,Sqrt(3)/2,y min ,y 1 ) 
if YHi < -1 then 
begin 

R := 1/i-Hi - l/Y) 

W := 2.5(y + (Sqrt(3)/2)i?)i? 2 

if W < Wends then return REJECT 

end 

7. F left := Sqrt(F /# ) 

if (Ymin < y> and Y Mt < Y max ) then Y Mt := Max(Yi, Yi eft ) 

else return NORESULT 
Y A := avgwt(yi cft ,y max ,.5) 
Z 2 := Sqrt(y max -y 2 ) 
Z 4 := -Sqrt(y max - y 4 ) 
if lOOOCi > 998 then 

begin 

T:= avgwt(y 1 ,y 2 ,l/16) 

A, := (T-y 1 )(33/16)+Integrate(Dx(-,^,F l ),T,y 2 ) 

A :=-(yieft-Yi)Sqrt(3) 

T:= avgwt(y cft ,y 4 ,l/16) 

A := A +Integrate(Dx(-,i7 ,F ),T,y 4 ) 

if A > Aj then return REJECT 

A G := A D +Integrate(Dx max (-, H , F ), Z 4 , Z 2 ) 

if A > Aj then return REJECT 

if leCi then return NORESULT 

end 

8. verify (Ci < Sqrt(3)/2 or y < y 2 ) 
T := Sqrt(y - Y min ) 

Z\:= Compare (Ci,Sqrt (3)/2, -T,T) 

Z 3 := Sqrt(y 2 -y min ) 

Ai := Integrate (Vx min (-, Hi, Fi), Z ± , Z 3 ) 

A :=Integrate (Dx(-, # , F ), Y 1: Y 4 ) 

if Aj < A then return REJECT 

A G := A D +Integrate(Dx max (-, H , F ), Z 4 , Z 2 ) 

if Aj / A then return REJECT 

9- WW* := Integrate(Dv min (-,iJj,Fi),Zi,Z 3 ) 

Wi := Wends + Wbase 

W :=Integrate (Dv(-, H , F Q ), Y 1 ,Y 4 ) + 

Integrate (Dv max (-, H , F Q ), Z 4 , Z 2 ) - W hase 
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if Wi / W then return REJECT 
else return NORESULT 

ymax 




Figure 17. Compared volumes. If we can show that V is not equal to 
Vbase + Kinds then we can discard the torus bubble as a potential equal 
volume minimizer. 

Theorem 5.3. If CheckRectangle(Ci,# ) returns REJECT without 
causing an IEEE exception, then there is no area- minimizing equal-volume 
torus bubble with cos 9\ € C\ and h Q S H Q . 

Proof. The procedure CheckRectangle is fed a pair of intervals containing 
a range of values for cos#i and h a , these intervals being denoted by C\ and 
H Q . We will call such a range of input values an input rectangle. If every pair 
of values (cos#i,/i ) 6 C\ x H a in an input rectangle can be ruled out as a 
possible minimizing torus bubble by some test, then the procedure returns the 
value REJECT. Otherwise it returns the value NORESULT (which results in 
a subdivision of the input rectangle into four smaller rectangles.) The tests 
applied to the rectangles of values for cos#i and h are based on the results 
developed in this paper, as we will explain below. At certain points we need 
to check whether some hypothesis required in our tests are satisfied. We can 
do so by using a verify statement. This statement will cause the program 
to terminate if it is not satisfied. The proof proceeds by interpreting and 
justifying each step of the algorithm. 

Step 1 rejects torus bubbles based on Proposition 4.21. G\ is a range of 
cosine values for the angle Q\, so lOOOCi > 996 corresponds to cos#i > .996, 
which implies that 9\ < 5.2° and the proposition applies. 

Step 2 first defines some values used in further steps. The values of F Q and 
Fi are calculated using the formulas in Proposition 4.16. It then examines 
the possible values of yi = sin 02- Corollary 4.10 implies —.5 < cos 62, and 
Lemma 4.8 implies cos #2 < -5. Since we assumed that 9\ < 62, we have 
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cos 02 < cos 9\ . Proposition 4.14 implies that — cos#i < cos 6*2- Thus the values 
for cos 62 arising from possible minimizing torus bubbles lie in the intersection 
of the intervals [—.5, .5] and [— cos 0\, cos 9\}. 
The force equation for Tj, 

fi(c) = (hi - 1)(1 - c 2 ) - Cv /3(l-c 2 ), 

has at most two roots by Lemma 4.15, one of which is c±. The other one is c 2 - 
Letting y = Vl — c 2 and rearranging gives 

fi = (hi - l)y 2 + cyV3 
(f t -(h t -l)y 2 ) 2 = 3y 2 (l-y 2 ) . 

This is a quadratic in y 2 . The sum of the two roots in y 2 is 

2(hj - l)h + 3 

3 + (hi-iy 

and the known root is 1 — Ci 2 , so subtraction gives the other root. Thus 
the interval T contains y 2 , the square of the second root. We next add the 
interval C2 2 , containing cos 2 6*2, to the interval T containing sin 2 62 to see if 
they contain values adding to one. If not, there is no minimizing torus bubble 
in the chosen input rectangle C\ x H Q , and this input rectangle is rejected. 

Step 3 calculates the possible values of X2 = cos 62 corresponding to y-values 
whose square is in the interval T. The only possible values for y2 are in [0, 1], so 
we intersect T with [0, 1] before taking the square root of T to get an interval 
y 2 which contains possible values of y^. Since 

fi = (hi - 1)(1 - C 2 2 ) - C 2V /3(1 - C 2 2 ) 

by Proposition 4.16, and 1 — C2 2 = y 2 , we can solve to get an expression for C2 
in terms of hi, fi and y2- 

C 2 = ((Hi - 1)Y 2 - F i /y 2 )/Sqrt(3) . 

Note that this expression is chosen to give us c 2 with the correct sign. We 
intersect the resulting interval with the interval C 2 computed in the previous 
step. If the resulting interval is empty, there is no minimizing torus bubble in 
the chosen input rectangle and it is rejected. 

Step 4 rejects torus bubbles based on Proposition 4.18. It uses the fact that 
cot (0i ) = ci/yi. 

Step 5 passes on a request to subdivide the rectangle if it is too wide. It 
does not reject anything, so there is nothing to justify. This step is included 
to speed the program's execution. 
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Step 6 calculates rough bounds to check if the volumes of the two regions 
can be equal. We calculate volume divided by tt in the algorithm to avoid the 
needless step of multiplying by the nonrepresentable constant it. The interval 
Wends = Kinds/ 7 '" is defined to be the volume of the solid of revolution between 
the spherical ends and the x-axis, divided by n. See Figure 17. This solid is 
wholly contained in the ball component. Y m \ n contains the y-coordinate of the 
local minimum of the Tj, based on solving the force equation as in Part 7 of 
Proposition 4.12. This minimum is realized in the torus bubble if and only if 
Q\ > 30°. Similarly Y mSuX contains the maximum height for r , as in Part 8 of 
Proposition 4.12. It is always achieved. Step 6 goes on to reject potential torus 
bubbles if a calculation shows that the torus component has less than half the 
total volume. The volume of the ball component is bounded below by the 
volume under the spherical caps, whose value is given by V^ n( js- The interval 

Y is defined to ensure that it contains the minimal y-coordinate of the torus 
component, which occurs either in Y m i n (if 9\ > 30°) or in Y±. An upper bound 

V for the torus component volume is calculated using the formula derived in 
Lemma 4.24. 

Step 7 tries to eliminate torus bubbles with small 9\, namely those with 
cos#i > .998, a condition which forces 9\ to be less than 3.7° and allows 
Lemma 4.26 to be applied. These small angles must be treated differently 
because the equations describing Delaunay curves are not well behaved near 
#1=0. Bounds on the x-displacement X2 — x\ are calculated using each of 
Tj and t . The interval A, contains an upper bound for the width of Tj and 
the interval A G contains a lower bound for the width of r . If the resulting 
intervals satisfy Aj < A G then Tj does not travel as far to the right as r and 
the input rectangle can be rejected. 

Step 7 begins by defining Y\ e f t to be an interval containing the y-value of 
the point where r G has a vertical tangency on the left as in Figures 18 and 19. 
This value is computed by taking the square root of the nonnegative part of 
F /H . Corollary 3.3 guarantees that this value is positive in a torus bubble, so 
we are justified in truncating negative values which may have appeared in the 
calculation before taking the square root. The outer curve is always a nodoid 
by Corollary 4.19, and if extended to a vertical tangency, the y-coordinate of 
that point is in the interval Y\ e ^. Replacing Yj e f t with Max(YJ e f t , Y\) assures 
that YJ e ft contains the y-coordinate of a vertical tangency if there is one, and 
the initial point otherwise. In particular it always contains the y-coordinate of 
the leftmost point. 

The value Y4 is used as an intermediate point for chopping up the inter- 
vals of integration for r G between its vertical and horizontal tangencies. See 
Figures 18, 19, 20. The limits of integration Z% and are defined for use 
when integrating near a maximum of r , as in Proposition 3.7. The square 
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roots taken in these calculations are justified in discarding the negative part of 
any interval, because for any particular torus bubble the quantity whose root 
is being taken is nonnegative. 




Figure 18. A double bubble configuration for < 6\ < 30. 



y, 




Figure 19. A double bubble configuration for 30 < 9\ < 60. 




Figure 20. A double bubble configuration for 60 < 9\ < 90. 

An upper bound for the x-displacement of Tj as it goes from y\ to y2 is 
given by Aj, and a lower bound for the x-displacement of r a as it goes from 
yi to Di is given by A Q . If Aj < A D , then a torus bubble cannot occur in the 
given input rectangle. 

We use the Intermediate Value Theorem to estimate Aj up to height y\ + 

— — — and A up to height y\ e ft- The estimates for Aj and A D near the x-axis 
16 

are justified by Proposition 3.4 in the case of r and by Lemma 4.26 in the 
case of Tj, as shown in Figure 21. 
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y 2 



upper bound for A; 



lower bound for A 

Figure 21. Estimating the widths of Tj from above and r G from below 
with 0\ small. 

For 0\ < 4°, Tj is strictly increasing, so that the x-displacement is equal to 
the integral of Dx from y\ to y 2 . This can be broken up into an integral from 
Y\ to T and an integral from T to Y 2 , where T is set to be one sixteenth of the 
way from Y\ to Y 2 . The values in T are all less than 1/2, so that Lemma 4.26 
applies and the slope of Tj below height T is smallest at its initial point (xi,yi). 
It follows that 

dx/dy < cot (30 - 0i) < cot 26° < 33/16 

on this interval and the integral from Y\ to T is bounded above by (33/16) 
(T - Yi). Then A» = (T - Yi)(33/16)+Integrate(Dx(-, Hi,Fi),T, Y 2 ) gives an 
upper bound for the width of Tj as it increases from Y\ to Y 2 . 

With r Q we first find an upper bound on the leftward displacement between 
j/i and yi c ft using a similar method. In this case the slope has absolute value 
larger than tan 30° = l/y/3 and (Y lcit - Yi)Sqrt(3) gives an upper bound for 
the leftward displacement to height yi e ft. An integral is calculated which gives 
a lower bound for the rightwards displacement of r G between yi e f t and y 2 , 
by integrating dx/dy from a point strictly above y\ e ft to y 2 - This integral is 
calculated in two pieces, the first to Y4 and the second from Y4 through the 
local maximum to Y 2 . A lower bound for A G is gotten by taking the upper 
bound for the leftwards displacement previously calculated and subtracting it 
from the sum of these integrals. If A > Aj the input rectangle is rejected. 

This step can reject torus bubbles even if the interval C\ contains the 
value 1, corresponding to 0\ = 0. If this step test fails to reject a range of 
potential torus bubbles, we check whether 1 is in the interval C\. If it is, then 
the next steps will not help, since they are ineffective for 0\ = 0. So if 1 is 
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in C\ we return NORESULT and send the current input rectangle back for 
subdivision. 

Step 8 calculates intervals for the x-displacement of Tj and r G in the case 
where the angle is larger than those treated in the previous step, namely 
cos 6*i < .996. In this step A Q is an interval containing the x-displacement 
of the outer curve and Aj an interval containing the x-displacement of Tj. In 
contrast to Step 7, they do not contain bounds for these displacements, but 
the actual widths of Tj and r D . Thus we can reject the input rectangle in this 
step if Aj and A D do not overlap. 

The interval A D is calculated by integrating Dx from Y\ to an intermediate 
I4 and then using Dx max to integrate from I4 across a maximum to Y^. If 
the width of A D is larger than 20 then the procedure returns NORESULT, 
calling for a subdivision of the input rectangle of intervals. This is done be- 
cause otherwise large amounts of computational time are wasted in inefficient 
calculations. 

The calculation of Aj is complicated somewhat because we need to deal with 
both the case where 6\ > 30° and the case where 9\ < 30°. Only in the former 
case, where C\ < Sqrt(3)/2, does Tj achieve an interior minimum. In that case, 
an integral along Tj can be calculated by integrating Dx m ; n from Z\ to Z%. If an 
interior minimum is not achieved, we can still integrate Dx m i n to calculate the 
x-displacement from Y\ to Y2, but we need to reverse the sign of Z\ as in the last 
statement of Proposition 3.7. Our integral formulas for x-displacement need to 
have that Y\ < Yi in this case, so we check with a verify statement that either 
this holds or there is an interior minimum. The Compare statement ensures 
the correct sign is given to Z\ unless C\ overlaps the interval Sqrt(3)/2. Z\ 
equals -Sqrt(Yi — F m ; n ) if Tj achieves an interior minimum and Sqrt(Yi — Y m i n ) 
otherwise. The appropriate sign is given to Z\ using a Compare statement. 
If C\ overlaps the interval Sqrt(3)/2 we can't tell which sign is correct, so the 
Compare statement returns an interval containing both Z\ and —Z\, and the 
integral then contains the x-displacements for both cases. If Aj and A G are 
disjoint, then the curves' final points cannot coincide, and the input rectangle 
is discarded. 

Step 9 calculates volumes divided by ir, as in Figure 17. We again calculate 
volume divided by tt in the algorithm to avoid an unnecessary multiplication 
of all volumes by tt. Vbase is the volume surrounded by Tj and Wbase = Vbase 1^. 
Adding Wbase to Wends gives Wi = Vi/ir). The calculation of V Q is slightly more 
complicated because the overhang on the left involves a subtraction, if there is 
an overhang. Such an overhang occurs when 9\ < 60°, as shown in Figures 18 
and 19. The formulas are set up so that volume is counted with a negative 
sign when r D is oriented to the left and positive sign when r is oriented to the 
right, as discussed in Proposition 3.6. This is just what is required to calculate 
the volume of the overhang. The volume Vbase is subtracted from the volume 
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inside T Q to get V , the volume of the torus component. Step 9 then rejects an 
input rectangle if the values W Q = V /ir and W{ = Vi/it are unequal. □ 

Now that we have tests that are able to reject certain ranges of torus bub- 
bles, we apply these tests to see whether they in fact reject all the possibilities. 
Since we don't know in advance how fat the ranges can be, we initially feed in a 
rectangle of values with < sin#i < 1,0 < h Q < 10. We recursively subdivide 
further as necessary. We monitor IEEE exception flags at this level, so that 
no computation is trusted if it raised an exception. 

Boolean procedure DivideAndCheckRectangle 
input : Y\, Hi 

Ci := Sqrt(l - Y?) 

call CheckRectangle (Ci, H Q ) 

if result is REJECT then return SUCCESS 

split Y\ in half, into Yi a ,Yib 

split H Q in half, into H oa , H Q b 

call DivideAndCheckRectangle ( Y± a , H oa ) 

call DivideAndCheckRectangle (Yi a , H b) 

call DivideAndCheckRectangle (Y15, H oa ) 

call DivideAndCheckRectangle (Yi&, H t>) 

return SUCCESS 

Main program 
begin 

clear exceptions 
begin 

Yx := [0, 1] 
H := [0, 10] 

Call DivideAndCheckRectangle (Yi, H ) 
end 

verify no exceptions raised 

Print ''All torus bubbles rejected. ' ' 

end 

Theorem 5.4. The algorithm described in Main, if run to completion 
without causing an exception, shows that no area minimizing torus bubble can 
enclose equal volumes. 

Proof. By Proposition 4.13, any minimizing torus bubble is determined by 
6*i and h Q , where 0\ is the angle of first arc and h Q is the mean curvature of T Q . 
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It suffices to consider < 9\ < 90° by Proposition 4.14, and < h Q < 10 by 
Proposition 4.3 and Proposition 4.25. 

Main calls DivideAndCheckRectangle with a rectangle containing the 
parameters of all possible torus bubbles. DivideAndCheckRectangle sub- 
divides this rectangle, depending on the success of its calculations. To avoid 
unnecessary use of trigonometric functions in the computation, the space of an- 
gles is parameterized by a variable y\ = sin#i. Thus the intervals Y\ cover the 
interval [0,1], in 1-1 correspondence with < Q\ < 90°. DivideAndCheck- 
Rectangle begins by using Y\ to compute C\, an interval corresponding to 
values of the cosine of an interval of 6\ angles. It then passes the C\ and H Q 
intervals to CheckRectangle where they are tested as previously described. 

Main can only finish if CheckRectangle rejects all rectangles. By Theo- 
rem 5.3, torus bubbles rejected by CheckRectangle are not minimizers. □ 

We implemented the algorithm described above in C++, though we could 
equally well have used any language which allows for control of rounding and 
follows the IEEE 754 standard, or any similar precise prescription for floating 
point arithmetic. We ran code implementing this algorithm on a variety of 
machines, various Wintel, Sun, HP, SGI, UNIX and Linux platforms, with 
identical results. The C++ code with embedded instructions is available at 
the web site [21]. The code takes about 10 seconds to run on a fast (1999) PC. 
The input rectangle was recursively subdivided up to ten times, resulting in a 
total of 15,016 rectangles examined. The calculation involved a total of 51,256 
integrals. 

From examining the output of the program, we can see the reasons for 
which various rectangle were rejected. Of the 15,016 rectangles examined, 
Step 1 rejected 31, Step 2 rejected 1140, Steps 3, 4 and 5 rejected none, Step 
6 rejected 3664, Step 7 rejected 541, Step 8 rejected 6691 and Step 9 rejected 
2949. In Figure 22 we indicated which steps were used to reject torus bubbles 
in which input rectangles. This chart offers some insight into why torus bubbles 
fail to be minimizers. However there is some arbitrariness to this, as doing the 
steps in a different order would have produced different data. 

In running the program no exceptions were raised, and the program returned 
the statement "All torus bubbles rejected." As a consequence we obtain a proof 
of our main result. 

Theorem 1. The unique surface of least area enclosing two equal volumes 
in M 3 is a double bubble. 

Proof Theorem 5.4 eliminates the possibility of a minimizing torus bubble. 
Theorem 4.1 implies that the only other possibility is a double bubble. 
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Figure 22. Chart showing which tests were used in each input rectangle. In 
some rectangles there were several subdivisions, and several tests were used. 
Steps 3, 4 and 5 were never used to reject input rectangles. 
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